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ABSTRACT 


Detailed  observations  of  wave  evolution  and  wave-eurrent  interaetion  in  tidal  inlets  and 
river  mouths  are  essentially  non-existent  owing  to  the  diffieulty  of  installing  and 
maintaining  fixed  instruments  in  this  harsh  environment.  This  work  develops  and 
explores  the  use  of  small,  free-drifting  buoys  to  eolleet  wave  and  eurrent  measurements 
in  eoastal  inlets.  The  instruments,  referred  to  as  Wave-Resolving  Drifters  (or  WRDs),  are 
small,  lightweight  and  inexpensive  enough  to  be  deployed  and  retrieved  in  large  numbers 
from  small  vessels.  To  study  wave  evolution  in  the  San  Franeiseo  Bight  and  the  Mouth  of 
the  Columbia  River,  30  WRDs  are  deployed  during  peak  ebb  tide  so  that  the  drifters  flow 
into  the  ineident  wave  field.  Wave  statisties  estimated  through  loeal  ensemble  averaging 
of  drifter  observations  and  ensemble-averaged  wave  speetra  are  used  to  deseribe  the  wave 
evolution  through  the  inlet  area.  The  observations  reveal  dramatie  spatial  variability  in 
the  wave  field  and  sometimes  doubling  of  the  loeal  wave  heights.  Comparisons  with 
numerieal  simulations  of  the  SWAN  (Simulating  Waves  Near  Shore)  model  and 
geometrie  opties  theory  (ray  diagrams)  show  the  distinet  effeets  of  refraetion  by  variable 
shoals  and  eurrents  on  the  wave  field  and  hint  at  nonlinear  instabilities  that  may  eause 
rogue  wave  development. 
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I.  INTRODUCTION 


Tidal  inlets  and  river  mouths  are  the  shipping  gateways  between  the  oeeans  and 
rivers  for  eommercial  vessels  and  the  U.S.  Navy.  Along  the  west  eoast  of  the  United 
States  two  of  these  major  gateways  are  the  Golden  Gate  in  the  San  Franeiseo  Bight  and 
the  Mouth  of  the  Columbia  River  where  the  interaetion  of  tidal  eurrents,  eomplieated 
shoal  systems  and  oeean  waves  ean  make  for  hazardous  navigation  eonditions. 

The  strait  under  the  Golden  Gate  Bridge  connects  the  San  Francisco  Bay  and  its 
complex  estuary  system,  which  is  the  largest  estuary  on  the  Pacific  Coast,  to  the  Pacific 
Ocean.  The  combined  economic  impact  of  the  Port  of  San  Francisco  and  the  Port  of 
Oakland  is  estimated  to  be  about  4  billion  dollars  each  year  (San  Francisco  Port 
Authority).  To  reach  these  two  ports,  more  than  4,000  ships  that  pass  through  the  Golden 
Gate  each  year  must  transit  along  a  shipping  channel  that  splits  the  bifurcated  shoal 
system  of  the  Four  Fathom  Bank  and  the  South  Shoal.  This  region,  known  as  the  San 
Francisco  Bight,  is  open  to  the  predominantly  west-northwest  swell  conditions  that  can 
pose  a  significant  hazard  to  navigation  when  they  steepen  in  response  to  the  shoal  system. 

San  Francisco  Bay  is  also  known  to  have  very  strong  ebb  tidal  currents  (about  2  m  s~’ ) 
that  interact  with  and  further  complicate  the  wave  field.  The  safest  navigation  path  is  to 
stay  in  the  channel,  but  that  too  can  be  dangerous  in  some  wave  conditions  (San 
Francisco  Port  Authority). 

An  estimated  3600  ocean-going  vessels  pass  through  the  Mouth  of  the  Columbia 
River  making  their  way  to  Portland,  OR  (Columbia  River  Bar  Pilots).  These  vessels  must 
pass  over  the  treacherous  region  known  as  the  graveyard  of  the  Pacific,  or  the  Bar,  as  the 
local  mariners  know  it.  When  westerly  approaching  ocean  waves  interact  with  the 
complicated  shoal  system  outside  the  mouth  of  the  Columbia  River  dangerous  wave 
conditions  can  occur.  The  natural  drainage  of  the  Columbia  River  couples  with  a  large 
semi-diurnal  tide  resulting  in  a  strong  ebb  tidal  jet  with  velocities  that  can  exceed  4  m  s”’ 
as  it  exits  the  shoal  system.  This  strong  current  enhances  the  already  rough  wave 
conditions  over  the  bar,  and  extreme  wave  conditions,  or  rogue  waves,  can  result  (Oregon 
State  Marine  Board). 
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This  dissertation  investigates  evolution  of  the  waves  in  the  approaches  to  the  San 
Francisco  Bay  and  the  Mouth  of  the  Columbia  River.  Traditional  measurement 
approaches  that  require  moored  instrumentation  are  helpful,  but  the  highly  variable  nature 
of  a  wave  field  interacting  with  a  complicated  shoal  system,  and  the  inherent  difficulty  of 
maintaining  fixed  instruments  in  a  shipping  channel  with  strong  currents  make  that 
approach  impractical.  Therefore,  a  fleet  of  50  wave-resolving  drifters  (WRD)  that 
measure  positions  and  horizontal  velocities  with  GPS  data-loggers  and  vertical 
accelerations  with  accelerometers  was  developed.  The  use  of  these  instruments  was 
explored  with  multiple  deployments  at  both  field  sites,  but  this  work  a  focuses  on  a  single 
deployment  from  each  site  that  best  highlights  the  new  application  of  this  technology  and 
identifies  the  physical  mechanisms  that  influence  the  wave  fields.  Chapter  II  outlines  the 
drifter  development,  deployment  and  analysis  technique  as  well  as  the  analysis  of  wave 
evolution  in  the  San  Francisco  Bight.  In  Chapter  III,  a  similar  analysis  is  applied  to  the 
Mouth  of  the  Columbia  River  where  extreme  waves  were  observed.  Finally,  Chapter  IV 
summarizes  the  findings  and  discusses  possible  future  work. 
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II.  DRIFTER  OBSERVATIONS  OF  WAVE  EVOLUTION 
IN  SAN  FRANCISCO  BIGHT 

A.  INTRODUCTION 

In  coastal  areas  characterized  by  strong  currents  and  energetic  waves,  sueh  as  in 
and  near  tidal  inlets,  in  situ  measurements  of  waves  and  currents  are  very  diffieult  to 
make.  In  particular,  the  use  of  standard  measurement  teehniques,  using  bottom-mounted 
or  moored  instruments,  is  eomplieated  due  to  e.g.  the  distortion  of  buoy  response  by 
eurrent  drag  on  mooring  lines  (Anetil  et  al.  1993;  Steele  1997;  de  Vries  et  al.  2003), 
burying  of  instruments  by  dynamic  seabed  morphology  (e.g.,  Barnard  et  al.  2006),  and 
loss  or  damage  of  instruments  by  ship  traffic  (Elias  et  al.  2012).  Moreover,  waves  and 
eurrents  generally  undergo  strong  spatial  variations  in  coastal  inlets  and  this  spatial 
evolution  and  inhomogeneity  is  difficult  to  capture  with  a  limited  number  of  fixed 
instruments.  As  a  consequenee,  synoptic  observations  of  wave-eurrent  interaction  and 
wave  evolution  in  such  regions  are  extremely  rare,  which  hampers  progress  in 
understanding  these  dynamics  in  natural  inlets. 

In  the  present  work,  we  explore  the  use  of  instrumented  drifters  to  resolve  wave- 
current  interactions  in  the  approaehes  to  and  inside  of  a  coastal  inlet.  Sueh  free-drifting 
instruments  can  be  readily  deployed  in  these  environments  and,  when  used  in  large 
numbers,  ean  provide  synoptic  information  on  the  variability  of  waves  and  currents  in  the 
area  of  interest.  The  observations  are  Lagrangian  (not  Eulerian),  and  in  regions  with 
strong  currents  and  spatially  variable  wave  fields,  recorded  time  series  can  be  highly  non¬ 
stationary  (even  if  the  wave  conditions  are  fairly  stationary  from  an  Eulerian  point  of 
view).  However,  by  deploying  a  large  number  of  drifters,  reliable  statistical  estimates  ean 
be  obtained  through  ensemble  averaging  of  statistieal  quantities  estimated  from  short 
time  series  of  nearby  drifter  measurements. 

In  the  past,  GPS  drifters  have  been  used  extensively  to  measure  oeean  currents 
and  surf  zone  dispersion  (MaeMahan  et  al.  2009),  but  recent  advances  in  GPS  teehnology 
have  made  it  possible  to  simultaneously  eollect  both  wave  and  current  measurements, 

even  with  inexpensive  GPS  reeeivers  (Herbers  et  al.  2012).  Although  the  absolute 
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position  accuracy  for  inexpensive  off-the-shelf  GPS  receivers  is  somewhere  between  2.5 
m  and  10  m,  these  errors  fluctuate  on  time  seales  much  longer  than  the  periods  of  wind 
waves  and  swell  and  thus  do  not  signifieantly  affect  wave  measurements  (Herbers  et  al. 
2012).  Moreover,  Doppler  veloeity  measurements,  whieh  are  also  available  on  most  GPS 
reeeivers,  are  far  less  sensitive  to  so  ealled  “jumps”  when  satellites  enter  and  leave  the 
useable  horizon.  As  a  eonsequenee,  using  the  Doppler  veloeities  rather  than  positions, 
generally  results  in  eleaner  signal  and  better  resolution  over  short  time  intervals,  whieh 
ean  be  particularly  important  in  the  wind  wave  band  (Thomson  2012). 

The  drifters  developed  in  this  study  utilize  the  Doppler  veloeity  measurements 
from  the  same  GPS  reeeivers  (Loeosys  GT-31)  as  those  used  by  Herbers  et  al.  (2012),  but 
the  sensor  paekage  is  augmented  with  an  aeeelerometer  to  resolve  vertieal  motions  and 
higher-frequeney  waves.  By  resolving  both  the  horizontal  and  vertical  wave  orbital 
motions,  we  remove  a  180-degree  ambiguity  in  direetional  estimates  based  on  horizontal 
GPS  measurements  alone  (see  Herbers  et  al.  2012).  The  proeedure  enables  the  aeeurate 
measurement  of  strongly  nonlinear  waves  (e.g.,  near  breaking  and  steep  waves  on 
opposing  eurrents). 

The  objeetive  of  the  present  work  is  twofold.  The  first  objeetive  is  to  validate  the 
improved  drifter  sensor  paekage  in  a  natural  (random)  wind  wave  field.  In  partieular  we 
test  the  ability  of  the  new  motion  sensor  paekage  that  eombines  the  GPS  and 
aeeelerometer  measurements  to  resolve  vertieal  motions  at  higher  frequeneies  than  was 
possible  with  the  drifters  used  in  Herbers  et  al.  (2012).  This  is  an  important  improvement, 
whieh  will  greatly  enhanee  the  potential  of  these  drifting  buoys  for  wave  observations. 
We  test  this  new  eapability  by  deploying  drifters  alongside  a  Datawell  Waverider  buoy  in 
relatively  homogenous  open  oeean  eonditions,  and  eompare  the  frequeney  and  directional 
spectra  of  the  different  buoys. 

The  seeond  objeetive  is  to  use  the  wave-resolving  drifters  (WRDs)  to  study  wave 
evolution  in  a  natural  inlet,  thus  exploring  their  use  in  regions  charaeterized  by  strong 
eurrents.  A  goal  is  to  develop  a  reliable  analysis  teehnique  that  ean  be  used  to  investigate 
the  physieal  proeesses  that  drive  wave  evolution  in  a  tidal  inlet.  To  this  end,  we  deployed 

large  numbers  of  WRD  during  ebb  eurrents  in  the  San  Franeiseo  Bight  (Figure  1)  to 
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observe  the  interaetions  of  wind  waves  and  swell  with  strong  opposing  eurrents.  The 
nonstationary  nature  of  the  data  records,  caused  by  the  rapid  transit  of  the  drifters  through 
the  channel  and  over  the  ebb  tidal  shoal,  complicates  the  analysis.  Here,  we  show  that 
robust  estimates  of  wave  spectra  and  bulk  statistics  can  be  obtained  from  these  records  by 
ensemble  averaging  short-time  estimates  over  many  drifters.  Additionally,  ray  diagrams 
and  simulations  from  a  SWAN  wave  model,  one-way  coupled  with  a  DelftSD  FLOW 
model,  are  used  to  identify  the  physical  mechanisms  responsible  for  the  changes  in  the 
wave  field  in  the  approaches  to  and  in  the  tidal  inlet  as  seen  in  the  drifter  observations. 

The  San  Francisco  Bight  field  experimental  site  and  drifter  sensor  package  design 
are  described  in  section  B.  The  wave-current  model  implementation  is  outlined  in  section 
C,  and  the  observed  wave  evolution  is  compared  to  model  predictions  in  section  D. 
Finally,  we  summarize  our  findings  in  section  E. 

B,  FIELD  SITE  AND  INSTRUMENTATION 

1.  Field  Site 

The  San  Francisco  Bay  is  a  large  tidal  estuary  that  is  connected  to  the  Pacific 
Ocean  by  a  narrow  inlet.  The  bay  side  (east)  of  the  inlet  is  spanned  by  the  Golden  Gate 
Bridge  at  the  narrowest  point,  and  the  west  side  of  the  inlet  opens  to  the  Pacific  Ocean.  A 
mixed  semi-diurnal  tidal  cycle  that  consists  of  a  28-day  spring  and  neap  cycle  (NOAA, 
2009)  results  in  strong  tidal  forcing  through  the  inlet  where  the  currents  can  exceed  2.5 
m  s”’  near  the  inlet  mouth  and  maintain  velocities  of  1  m  s”’  to  the  ebb  tidal  shoal 
approximately  8  km  offshore  (Barnard  et  al.  2006;  Elias  et  al.  2012;  and  Barnard  et  al. 
2012).  Waves  approaching  the  San  Erancisco  Bight  from  a  predominantly  west- 
northwest  direction  interact  with  a  complex  shoal  area  bathymetry  and  a  strong  ebb  tidal 
jet  (Elias  and  Hansen  2012  and  references  therein). 

2.  Wave  Resolving  Drifters 

The  drifters  developed  and  used  in  this  study  were  designed  as  lightweight 
instrument  packages  that  can  be  deployed  from  small  vessels  and  accurately  resolve  both 
ocean  surface  waves  and  currents.  The  drifter  itself  consists  of  a  30  cm  diameter,  hard 
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shell,  closed-cell  foam  core  buoy  (Jim-Buoy  model  4400-RF)  with  a  1 .25  cm  wide  rubber 
band  attached  near  the  waterline  to  mitigate  the  immersion  buoyancy  resonance  that 
occurs  at  1.11  Hz  for  this  drifter.  A  ballast  weight  (10.4  kg)  is  suspended  from  a  shackle 
to  submerge  the  buoy  approximately  to  the  centerline  (see  Figure  1). 

The  wave-current  sensor  package  attached  to  the  buoy  consists  of  a  GPS  receiver 
(Locosys  GT-31),  and  a  three-axis  accelerometer  package  (Gulf  Coast  Data  Concepts  X6- 
2).  The  latter  is  tethered  between  the  buoy  and  the  ballast  chain  (Figure  1)  to  maintain  a 
nearly  vertical  orientation  of  the  accelerometer,  irrespective  of  the  buoy  orientation.  This 
design  suppresses  the  sensitivity  to  roll  and  pitch  motions  of  the  buoy  vertical 
acceleration  measurements  that  can  result  from  steeper,  high-frequency  waves. 

To  allow  real-time  tracking  of  the  drifter’s  position,  which  is  essential  as  a 
recovery  aid  when  deploying  large  numbers  of  drifters  close  to  shore  in  strong  currents, 
the  drifter  is  equipped  with  a  Garmin  DC  40  GPS  transmitter. 

a.  GPS  Position  and  Velocity  Sensor 

The  Locosys  GT-31  GPS  receiver,  based  on  the  SiRF  111  chipset,  collects  data  at 
IHz  and  accurately  resolves  the  horizontal  components  of  the  wave  orbital  displacements 
using  differential  position  data  (see  Herbers  et  al.  2012).  Additionally,  estimates  of 
horizontal  velocity  can  be  extracted  from  the  Doppler  shifts  in  the  raw  LI  carrier  phase 
signals.  The  resolution  of  the  Doppler  velocity  data  is  an  order  of  magnitude  higher  than 
the  differential  positioning  data,  which  particularly  improves  the  accuracy  of  high- 
frequency  wind-wave  measurements  (Thomson  2012).  In  these  experiments  both 
horizontal  Doppler  velocity  and  differential  positions  were  recorded  and  analyzed. 

Herbers  et  al  (2012)  show  that  the  GT-31  sensors  yield  accurate  and  reliable  wave 
energy  and  direction  spectra.  However,  vertical  motions  cannot  be  accurately  resolved 
with  these  inexpensive  GPS  receivers.  Although  the  horizontal  motions  alone  can  be  used 
to  estimate  wave  statistics  (as  shown  by  Herbers  et  al.  2012),  this  requires  the  use  of 
linear  wave  theory,  which  is  less  reliable  in  very  steep  waves,  and  wave  direction 
estimates  based  on  horizontal  motions  alone  suffer  from  a  180-degree  ambiguity  (Herbers 
et  al.  2012). 
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b.  Accelerometer 

To  augment  the  GPS  sensor  and  resolve  the  vertieal  motions,  a  three-axis 
aecelerometer  (X6-2  MEMS)  manufaetured  by  Gulf  Coast  Data  Coneepts,  EEC  was 
added  to  the  sensor  paekage.  The  aeeelerometer  measures  vertieal  aecelerations  in  the 

range  of  +/-  2  g,  where  g  is  the  gravitational  aoeeleration  (9.81  m  s“^) .  The  sampling 

_2 

rate  is  10  Hz,  and  aeeeleration  fluctuations  as  small  as  0.001  ms  are  resolved,  which 
make  this  accelerometer  well  suited  to  the  measurement  of  surface-wave  accelerations 
that  are  typically  a  small  fraction  of  gravity  (g).  Since  the  X6-2  is  a  stand-alone 
accelerometer  and  not  part  of  an  inertial  measurement  unit  (IMU),  it  does  not  resolve  the 
body  orientation  so  that  the  vertical  accelerations  may  contain  errors  due  to  changes  of 
the  pitch  and  roll  of  the  sensor  during  the  motion.  However,  these  errors  are  likely  small 
as  the  sensor  was  tethered  between  the  buoy  and  the  ballast  chain,  thus  suppressing  the 
pitch  and  roll  motions  induced  by  the  sea  surface  slope. 

The  accelerometer  wave  measurements  are  noisy  at  low  swell  frequencies  where 
the  signal-to-noise  ratio  suffers  from  small  accelerations  (low  steepness  of  these  longer 
waves)  and  contamination  by  nonlinear  buoy  motions  (Appendix  A). 

c.  Integrated  Sensor  Package 

To  account  for  the  different  sample  rates,  the  accelerometer  and  GPS  time  series 
were  interpolated  onto  a  common  10  Hz  time  base  (see  McIntyre  2013  for  more  details). 
The  accelerometer  is  capable  of  resolving  high-frequency  motions,  but  its  accuracy  is 
degraded  at  lower  frequencies  (<  0.15  Hz).  On  the  other  hand,  the  GPS  receiver  is  most 
accurate  in  the  swell  frequency  band  (-0.05  -  0.15  Hz)  but  becomes  increasingly  noisy  at 
higher  frequencies  (Herbers  et  al.  2012).  Therefore,  to  obtain  an  accurate  measurement  of 
the  wave  statistics  and  the  wave  spectrum  across  the  complete  frequency  range,  covering 
both  the  swell  and  wind  wave  bands,  we  estimate  a  composite  spectrum  that  combines 
the  GPS  (lower-frequency  range)  and  accelerometer  (higher-frequency  range). 
Specifically,  we  use  the  GPS  observations  in  the  range  0.05  Hz  -  0.15  Hz  and  vertical 
accelerations  (from  accelerometer)  in  the  range  0.25  Hz  -  1  Hz.  In  the  intermediate  range 
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(0.15  Hz  -  0.25  Hz)  the  spectrum  is  determined  from  a  weighted  average  of  the  two 
sensors,  where  the  GPS  (accelerometer)  weighting  was  varied  linearly  from  1(0)  at  0.15 
Hz  to  0(1)  at  0.25  Hz.  The  swell  peak  is  well  resolved  in  the  combined  sensor  spectrum 
(Figure  2),  in  good  agreement  with  the  Datawell  estimate,  and  the  high-frequency  tail 
closely  follows  the  f~^  (Phillips  1984)  high-frequency  slope  out  to  the  1.11  Hz 
maximum  frequency  that  can  be  resolved  by  the  Wave  Resolving  Drifter  (WRD)  buoy.  In 
Figure  2  the  Datawell  spectrum  is  shown  only  out  to  0.32  Hz  because  of  its  lower 
sampling  rate  (1.28  Hz)  and  immersion  resonance  frequency  (0.7  Hz).  Additionally, 
estimates  from  the  WRD  are  within  3%  of  the  Datawell  buoy  estimates,  which  was  used 
for  ground  truth  in  the  validation.  A  detailed  discussion  of  the  methodology  and 
validation  is  provided  in  Appendices  A  and  B. 


C.  WAVE-CURRENT  MODEL  IMPLEMENTATION 


To  support  the  analysis  of  the  observations  collected  in  the  San  Francisco  Bight,  a 
region  characterized  by  complicated  bathymetry  and  strong  currents,  we  implemented  the 
third-generation  wave  model  SWAN,  one-way  coupled  to  the  Delft3D  circulation  model 
for  this  area  (see  Elias  and  Hansen,  2012).  The  wave  model  SWAN  is  based  on  the  action 
balance  equation  (Booij  et  al.  1999  and  Ris  et  al.  1999),  which  can  be  expressed  as 


dN  , 


3(7 


(7 


(1) 


where  the  source  terms  can  be  written  as 


+  Snl4  +  ^bf  +  ^brk 


(2) 


The  left  side  of  equation  (1)  represents  the  conservation  of  the  action  density 
E 

N(G,d;  x,y,t)  =  — ,  where  E  is  energy,  o  is  relative  frequency  (moving  reference 

(7 

frame),  6  is  wave  direction,  and  x,y,  and  t  represent  the  horizontal  coordinates  and 
time,  respectively.  The  current  field  is  represented  by  U ,  is  the  (linear)  group  speed, 
and  Cg  and  represent  transport  rates  in  directional  and  frequency  space,  respectively 
(see  Booij  et  ah,  1999). 
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The  source  term  (the  right  side  of  equation  (1))  settings  are  as  follows.  The  wind 
input,  S-^,  was  set  to  the  default  3’^'*  generation  Cavaleri  and  Malanotte-Rizzoli  (1981) 
parameterization,  the  whitecapping  term  was  set  to  Van  der  Westhuysen  et  al.  (2007), 
nonlinear  quadruplet  interactions,  S^, ,  were  computed  with  the  standard  Discrete 
Interaction  Approximation  (Hasselmann  et  al,  1985),  bottom  friction  S^j-  was  set  to 

Hasselmann  et  al.  (1973)  using  the  swell  setting,  and  the  Battjes  and  Janssen  (1978) 
formulation  was  used  for  depth-induced  wave  breaking  . 

The  SWAN  simulations  were  done  in  nonstationary  mode  and  finite  differences 
were  solved  on  an  unstructured  grid  with  minimum  resolution  offshore  and  maximum 
resolution  in  the  region  of  interest.  To  accurately  resolve  refraction  effects  of  the  swell 
band  over  the  complicated  medium  variations  in  this  area,  we  applied  a  2-degree 
resolution  for  the  wave  directional  space  in  our  model  implementation.  The  offshore 
boundary  conditions  were  obtained  from  the  spectral  estimates  from  wave  measurements 
made  by  the  Point  Reyes  NOAA  buoy  #46214  (Figure  3). 

The  wind  forcing  was  separated  into  two  regions.  The  wind  measurement  time 
series  from  the  NOAA  San  Francisco  Buoy  #46026  (located  33  km  offshore)  was  applied 
uniformly  in  the  domain  seaward  of  the  mouth  of  the  inlet,  and  the  time  series  from  the 
NOAA  station  #FTPC1  located  just  inside  the  inlet  on  the  south  of  the  channel  was  used 
for  the  inlet  (small  box  in  Figure  3).  The  wind  measurements  from  these  locations  were 
found  to  be  reasonably  representative  (within  3  m  s”’  on  average)  of  the  measurements 
made  by  the  shipboard  anemometer  while  transiting  these  regions  at  coincident  times. 

The  currents  were  simulated  using  the  Delft3D  flow  model  using  the 

implementation  for  the  San  Francisco  Bay  developed  by  Elias  and  Hansen  (2012),  which 
was  developed  to  investigate  sediment  transport  in  regions  around  the  San  Francisco  Bay 
and  Bight  (see  Elias  and  Hansen  2012;  Barnard  et  al.  2012).  The  circulation  model 
consists  of  six  decomposed  (2-way  coupled)  domains  that  allow  for  varying  resolution, 
and  results  in  a  spatial  resolution  of  about  50  m  in  our  region  of  interest  (the  San 
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Francisco  Bight).  Due  to  the  extensive  seale  of  the  full  domain,  this  model  was  run  in 
depth-averaged  mode  to  gain  eomputational  effieieney  (Elias  and  Hansen  2012). 

The  SWAN  model  domain  used  in  this  study  extends  50  km  north  and  south  of 
the  Golden  Gate  and  50  km  offshore  to  the  Point  Reyes  buoy  (see  Figure  3).  A 
eomparison  of  the  speetra  and  signifieant  wave  height  {H ^  as  defined  in  Appendix  A) 
observed  at  the  San  Franeiseo  NOAA  Buoy  just  outside  our  region  of  interest  shows 
reasonably  good  agreement  with  the  SWAN  simulation  (Figure  4).  A  synoptie  view  of 
the  eurrent  field  and  predieted  signifieant  wave  heights  at  peak  ebb  flow  (Figure  5)  shows 
the  expeeted  strong  variability  of  waves  with  amplifieation  over  the  bar  and  in  the  eurrent 
jet. 

To  isolate  the  effeets  of  different  physieal  meehanisms  and  foreing  on  the 
observed  wave  evolution,  various  SWAN  simulations  were  eondueted  with  alternately 
wind  foreing,  eurrent  or  bottom  frietion  turned  off. 

D,  RESULTS 

1.  Drifter  Observations 

To  study  wave  evolution  in  the  San  Franeiseo  Bight  during  an  ebb  tide,  we 
deployed  an  array  of  30  WRDs  near  the  Golden  Gate  on  April  27,  2012  (see  Figure  1). 
The  WRDs  were  plaeed  near  the  eenter  of  the  ebb  tidal  jet,  so  that  (with  few  exeeptions) 
they  were  earned  aeross  the  ebb  tidal  shoal  (labeled  SF  Bar  in  Figures  1  and  3)  through 
the  opposing  wave  field.  The  large  array  of  drifters  allowed  for  ensemble  averaging  to 
estimate  wave  speetra  and  bulk  statisties  from  the  non-stationary  time  series  eolleeted  by 
the  drifters  as  they  traversed  the  spatially  inhomogeneous  wave  field  aeross  the  inlet. 

To  eapture  the  ebb  tidal  eyele,  the  deployments  started  about  1.5  hours  before 
predieted  peak  ebb  eurrent  (1.88  m  s”’),  near  the  seaward  entranee  to  the  San  Franeiseo 
Bay  inlet.  They  were  deployed  in  elusters  of  five  under  the  Golden  Gate  Bridge  with  one 
minute  spaeing  between  eaeh  drifter  in  a  eluster  and  six  minutes  between  eaeh  eluster. 
The  last  WRD  was  deployed  near  the  time  of  peak  ebb  eurrent,  whieh  ensured  that  all  the 
WRD  would  depart  the  mouth  of  the  inlet  as  the  ebb  tide  began  to  slaeken. 
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Since  the  drifters  were  earried  by  the  fast  flowing  ebb  eurrent  (surface  flows 
approaehing  2  ms”’)  through  a  spatially  inhomogeneous  wave  field,  the  time  series 
reeorded  by  the  individual  drifters  are  non-stationary.  To  enable  the  use  of  standard  time 
series  analysis  techniques  (assuming  stationarity),  we  divide  the  observations  in  relatively 
short  segments  (4096  samples  at  10  Hz  or  roughly  6.8  minutes),  and  to  improve  statistieal 
reliability,  we  average  speetral  estimates  over  ensembles  of  drifters  that  pass  through  the 
same  area. 

The  nearly  uniform  lateral  expansion  of  the  drifters  as  they  exit  the  channel  (see 
Figure  6a),  allows  us  to  estimate  wave  statisties  through  spatial  averaging  onto  a  regular 
500  m  X  500  m  grid  and  thus  obtain  a  synoptie  view  of  the  wave  variability  in  the  area  of 
interest.  Figure  6a  shows  the  ensemble-averaged  H ^  for  eaeh  grid  eell  plotted  on  the  data 
eentroid  of  the  drifter  ensemble  loeations.  The  number  of  drifter  observations  per  grid 
eell  varies  aeross  the  spatial  domain,  and  the  resulting  ensemble  sizes  vary  from  a  single 
drifter  estimate  at  certain  points  on  the  outer  edges  of  the  drift  tracks,  to  up  to  40  in  the 
well-traveled  areas  in  the  middle  of  the  jet. 

Over  the  San  Franeiseo  Bar  (SF  Bar),  a  dramatie  inerease  in  wave  height  (loeally 
up  to  3  m)  is  followed  by  a  fairly  abrupt  deerease  to  about  1.6  m  immediately  east  of  the 
SF  Bar.  Approaching  the  inlet,  the  wave  heights  again  inerease  as  the  waves  eneounter 
the  strong  tidal  jet,  followed  by  a  dramatie  deerease  inside  the  inlet.  A  similar  wave 
height  pattern,  albeit  with  smaller  with  more  gradual  variations,  is  defined  by  the  SWAN 
simulation  (Figure  6b). 

Additionally,  ensemble-average  speetra  were  eomputed  in  five  regions  where  the 
wave  charaeteristics  in  the  WRD  observations  are  similar.  The  regions  (see  Figure  6a) 
were  defined  as  follows:  Region  A  (offshore,  ensemble  average  of  489  spectra).  Region 
B  (SF  Bar,  273  speetra).  Region  C  (east  of  the  SF  Bar,  157  speetra).  Region  D  (approaeh 
to  the  inlet  entranee,  270  speetra),  and  Region  E  (inside  inlet,  172  speetra). 

In  the  discussion  that  follows,  we  will  focus  the  analysis  on  speetral  bands, 
defined  here  as  the  swell  band  (0.05  Hz  -  0.15  Hz),  the  wind  sea  band  (0.15  Hz  -  0.3 
Hz),  and  the  speetral  tail  (0.3  Hz  -  1.0  Hz).  We  will  eompare  observed  variability  of 
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these  bands  to  SWAN  simulations,  while  varying  the  various  foreing  terms  in  the  model 
simulations  (to  isolate  the  effeet  of  eaeh  foreing  term),  eombined  with  ray  eomputations, 
to  identify  the  dominant  physieal  meehanisms  that  eause  the  observed  variability  in  this 
region. 

The  discussion  in  each  section  will  progress  from  offshore  into  the  inlet,  thus 
following  the  wave  field  (as  opposed  to  the  drifters)  as  it  evolves  through  the  inlet. 

2,  Swell  Band 

The  WRD  estimates  ofH  ^  (dominated  by  '-ll  s  swell)  increase  progressing  from 
offshore  to  over  the  SF  Bar  (Figure  6a),  and  a  similar  increase  is  apparent  in  the  default 
SWAN  Sim  #1  (Figure  6b).  There  is  a  sharp  amplification  in  the  swell  band  of  the 
observed  and  predicted  spectral  estimates  (Figure  7a,  b),  indicating  that  the  wave  height 
variability  in  this  region  is  predominantly  related  to  swell  focusing  on  the  bar.  The  effects 
of  shoaling  due  to  the  decrease  in  group  speed  c^only  account  for  about  10%  of  the 

amplification  in  wave  height.  The  main  cause  for  the  wave  height  increase  is  the  strong 
refractive  focusing  that  is  evident  in  the  ray  trajectories  (Figure  8).  A  comparison  of  the 
rays  computed  with  currents  (Figure  8a)  and  without  currents  (Figure  8b)  shows  little 
change  and  suggests  that  the  focusing  is  predominantly  in  response  to  the  refractive  effect 
of  variable  water  depth  over  the  bar.  In  SWAN  Sim#4  (Figure  9c)  only  marginally  less 
amplification  occurs  over  the  bar,  which  also  suggests  that  the  amplification  is  due  to 
bottom  refraction. 

Continuing  east  across  the  bar,  the  observed  abrupt  decrease  in  (see  Figure  6a) 
is  also  evident  in  all  SWAN  simulations.  However,  this  decrease  is  most  pronounced  in 
Sim  #2  where  the  disabling  of  wind  forcing  accentuates  the  swell  wave  variability.  The 
abrupt  decrease  of  nearly  Im  in  wave  height  (Figure  6a)  and  the  similar  drop  in  swell 
spectral  levels  seen  in  both  the  WRD  and  SWAN  results  (Figures  7a  and  7b)  cannot  be 
explained  by  bottom  friction  alone  as  the  waves  transit  the  bar,  which  is  confirmed  by 
Sim  #3  (Figure  9b)  where  bottom  friction  is  disabled  yet  a  similar  reduction  in  wave 
height  occurs.  The  orientation  of  the  north  bar  and  the  SF  Bar  relative  to  the  prevailing 

swell  arriving  from  the  WNW  results  in  strong  refractive  effects.  The  dramatic  decrease 
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of  energy  in  the  swell  band  in  this  region  is  assoeiated  with  the  divergenee  of  the  ray 
trajeetories  in  the  lee  of  the  bar  (Figure  8a).  These  strong  refraetion  effeets  are  also 
evident  in  the  directional  spread  estimates  from  the  WRD  (Figure  10)  that  increase  from 
less  than  20  degrees  seaward  of  the  bar  to  greater  than  30  degrees  behind  the  bar.  A 
comparison  of  the  rays  with  currents  (Figure  8a)  and  without  currents  (Figure  8b)  as  well 
as  SWAN  Sim  #1  (Figure  6b  with  currents)  and  Sim  #4  (Figure  9c  without  currents), 
shows  that  the  effects  of  the  current  field  on  the  swell  is  negligible.  This  shows  that  the 
large  difference  in  wave  height  observed  between  the  bar  and  to  the  region  just  to  the  east 
of  it,  is  principally  due  to  refractive  focusing  by  the  bathymetry. 

Progressing  east,  a  slight  amplification  is  observed  in  Region  D  (Figure  6a), 
which  -  to  a  lesser  degree  -  is  also  evident  in  the  SWAN  simulations  (e.g.  Sim  #2  without 
wind.  Figure  9a).  However,  the  limited  increase  of  spectral  energy  levels  in  the  swell 
band  (Figure  7a)  does  not  account  for  the  overall  increase  in  as  evident  in  Figure  6a, 
which  suggests  that  changes  in  the  other  frequency  bands  are  important  here  as  well 
(which  is  discussed  in  the  following). 

When  progressing  farther  into  the  channel,  the  observed  and  predicted  H ^  again 
decrease  abruptly  (Figure  6a).  The  decay  is  especially  strong  in  the  swell  band  (Figure  7a, 
b)  and  appears  to  be  caused  by  refraction  due  to  varying  depth  (note  the  ray  divergence 
over  the  deeper  channel  in  Figure  8a,  b),  although  diffractive  effects  of  the  channel  (Melo 
and  Guza  1991),  which  are  not  accounted  for  in  SWAN,  may  also  be  important.  Although 
the  refraction  on  the  jet-like  ebb  current  will  likely  oppose  somewhat  the  refractive 
divergence  induced  by  the  bathymetry  in  this  region,  the  (longer)  swell  waves  are  more 
sensitive  to  variations  in  the  bathymetry  than  the  current  field. 

3,  Wind-wave  Band 

In  this  section  we  focus  on  the  evolution  of  less  energetic,  higher-frequency 
waves  using  the  average  spectra  in  regions  A-E.  Offshore  of  the  SF  Bar,  the  wind  wave 
band  (Figure  7a)  is  characterized  by  a  slight  elevation  in  energy  near  0.2  Hz  that  extends 
to  0.3  Hz  before  it  rolls  off  into  the  f~*  spectral  tail  (similar  to  the  spectrum  from 

homogenous  open  ocean  conditions  in  Figure  2).  The  currents  in  the  offshore  region  are 
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less  than  0.5  ms’,  but  increase  to  1.4  m  s  ’  over  the  SF  Bar.  A  noticeable  increase  in 
energy  in  the  wind  wave  band  is  apparent  in  the  spectra  (Figure  8a)  which  results  in  an 

Hk„ 

increase  in  the  non-dimensional  mean  steepness  parameter,  e  =  — ,  (Banner  2000) 

from  0.04  offshore  to  0.07  over  the  bar  (where  H ^  in  the  frequency  range  from  0.15  Hz 
to  1.00  Hz  is  considered  and^^  is  the  peak  wave  number  that  follows  from  the  linear 

dispersion  relation  at  the  relative  frequency  peak).  The  minimum  depth  on  the  shoal  is  8 
m  and  the  ray  trajectories  indicate  that  the  depth  variability  does  not  much  affect  the 
propagation  of  the  (relatively  short)  dominant  wind  waves.  However,  current  refraction 
(Figure  10a)  causes  a  slight  convergence  on  the  northeast  portion  of  the  SF  Bar. 

Contrary  to  what  was  observed  in  the  swell-band  (where  a  strong  decrease  was 
seen),  the  region  east  of  the  bar  exhibits  a  marked  increase  in  energy,  in  the  wind-wave 

band  (Figure  7a).  The  surface  velocity  in  this  region  increases  to  a  mean  of  1.88  m  s”’ 
resulting  in  further  steepening  (mean  steepness  parameter  of  e  =  .ll)  as  the  waves 
steepen  on  the  current  (Bascheck  et  al.  2006)  and  ray  trajectories  including  the  currents 
(Figure  11a)  show  pronounced  focusing  in  this  region  of  amplified  wind  wave  energy. 
Rays  without  currents  (Figure  11b)  show  virtually  no  focusing  in  this  region  (or 
anywhere  else  for  that  matter).  These  results  indicate  that  the  waves  in  this  band  are 
refracting  onto  the  strengthening  ebb  tidal  jet,  which  explains  the  increase  in  energy,  and 
steepening  ( e  =  0.1 1 )  observed  in  this  region. 

Approaching  the  entrance  to  the  inlet,  the  energy  levels  in  the  wind-wave  band 
continue  to  increase  (Figure  7a),  and  despite  the  marked  decrease  in  swell  energy  there  is 
an  increase  in  overall  in  this  region,  as  wave  conditions  change  from  swell- 
dominated  to  wind-wave  dominated.  The  ebb  tidal  jet  is  very  focused  in  this  region  (see 
Figure  5)  and  the  surface  currents  reach  a  maximum  mean  velocity  (1.89  ms”’). 
Additionally  the  north  shoal  (where  the  shallowest  depths  are  around  4  m)  refracts  the 
wind  sea  waves  into  the  channel,  where  they  are  superposed  by  waves  directly  incident 
onto  the  inlet  mouth  and  enhanced  by  the  current  jet  (compare  Figure  11a  and  11b).  The 
mean  steepness  parameter  in  this  region  reaches  0.15,  which  indicates  a  high  probability 
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of  breaking  (Banner  2000).  Current-induced  refraction  on  the  ebb  tidal  jet  and  steepening 
on  the  opposing  current,  results  in  a  region  of  enhanced  directional  spreading  and 
steepness  in  this  band. 

Finally,  entering  further  into  the  inlet,  spectral  energy  levels  drop  dramatically,  in 
part  attributed  to  refraction  as  ray  trajectories  for  both  swell  and  wind-wave  frequencies 
show  a  strong  divergence  in  the  channel  (Figures  8a,  11a).  Breaking  dissipation  of  high- 
frequency  waves,  steepening  on  the  opposing  current,  may  also  contribute  to  the  observed 
decay.  Additionally,  the  surrounding  terrain  provides  considerable  sheltering  in  this 
region,  resulting  in  a  dramatic  decrease  in  wind  forcing  conditions  within  the  inlet  itself 
(ship  observations  were  below  4  m  s~'  in  the  inlet),  virtually  eliminating  any 
regeneration  of  the  wind  waves  on  the  strong  current. 

4,  High-frequency  Tail 

Moving  from  offshore  to  the  SF  Bar  (region  B  in  Figure  8),  the  high-frequency 
tail  shows  a  slight  decrease  in  the  slope  at  frequencies  above  the  local  (linear)  blocking 

g 

frequency  (0.54  Hz,  estimated  from  f  = - ).  This  decrease  is  not  observed  in  the 

AnU 

SWAN  simulations  where  the  parameterized  f~^  slope  persists. 

Progressing  inside  the  bar  (Region  C),  the  slope  of  the  high-frequency  tail  begins 
to  deviate  from  f~^  at  lower  frequencies  presumably  in  response  to  the  decreased 
blocking  frequency  (0.42  Hz)  associated  with  the  stronger  current  in  this  region.  Decrease 
in  high-frequency  energy  levels  in  this  region  is  indicative  of  wave  breaking  dissipation. 
This  is  consistent  with  an  enhanced  mean  steepness  parameter,  e  =  0.11,  which  is 
significantly  larger  than  the  critical  value  e  =  0.05  suggested  by  Banner  (2000). 

Approaching  the  mouth  of  the  inlet  (region  D  in  Figure  8),  a  further  decrease  is 
observed  in  slope  of  the  high-frequency  tail  further  decreases  that  is  absent  in  the  SWAN 
results,  which  suggests  enhanced  dissipation  of  high-frequency  waves  is  occurring  as 
they  steepen  in  response  to  the  current. 
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Finally,  within  the  inlet  (region  E),  the  deerease  in  the  high-frequency  spectral 
levels,  visible  in  both  the  WRD  and  SWAN  spectral  estimates  (Figures  7a  and  7b)  may 
be  caused  by  blocking  and  dissipation  on  the  opposing  current,  preventing  these  short 
waves  to  enter  this  region.  This  combined  with  much-reduced  wind  velocities  due  to 
sheltering  effects,  eliminates  the  re-generation  of  high-frequency  waves  in  this  region. 

E.  SUMMARY  AND  CONCLUSIONS 

Strong  tidal  currents,  energetic  waves,  and  dynamic  seabed  morphology  can  make 
it  very  difficult  to  collect  observations  of  waves  and  currents  in  coastal  inlets  using 
conventional,  bottom-mounted  instruments.  As  a  consequence,  such  observations  are 
exceedingly  rare.  Moreover,  it  is  difficult  to  capture  the  regional  variability  and 
inhomogeneity  typically  seen  in  coastal  inlets  using  an  array  of  fixed  instruments.  In  this 
work  we  explore  the  use  of  free-floating  Lagrangian  drifters  as  an  alternative  in  situ  sensing 
strategy  for  measuring  waves  and  surface  currents  in  tidal  inlets.  By  treating  the  drifter 
observations  as  independent  realizations  and  ensemble  average  across  them  to  obtain 
accurate  statistical  estimates,  we  are  able  to  study  spectral  evolution  across  the  inlet  and 
identify  the  various  physical  mechanisms  that  most  influence  the  spatial  variability  of  the 
wave  field  across  various  frequency  ranges  (swell,  wind  sea,  high-frequency  tail). 

The  wave-resolving  drifters  (WRD)  combine  a  GPS  receiver  that  records  horizontal 
displacements  and  Doppler  velocities,  with  a  three-axis  accelerometer  to  resolve  the 
vertical  motions.  The  accelerometer  extends  the  WRDs  ability  to  resolve  the  high- 
frequency  tail  of  the  wind  wave  spectrum,  and  removes  the  180-degree  ambiguity  of  GPS 
wave  direction  estimates  (see  Herbers  et  al.  2012).  Comparison  to  a  conventional  Dataweb 
DWR-G  Waverider  buoy  shows  good  agreement,  with  slightly  elevated  energy  levels  away 
from  the  peak  resulting  in  a  positive  bias  of  approximately  3%  in  significant  wave  height 
estimates  (Appendix  A  and  B).  The  mean  directions  agree  to  within  1  degree,  but  the 
WRD  directional  spread  is  approximately  10  degrees  greater  than  the  Dataweb  estimates 
(Figure  2).  This  validation  confirms  that  the  WRD  developed  in  this  study  can  provide 
reliable  routine  heave-pitch-rob  type  wave  measurements  and  perform  at  a  similar  level  as 
the  more  expensive  and  larger  3-D  GPS  drifters  evaluated  by  Herbers  et  al.  (2012). 
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The  drifters  developed  here  are  relatively  low  eost  (<  $1000).  Henee,  they  can  be 
used  in  large  numbers,  which  allows  for  improved  estimates  of  wave  spectra  and  bulk 
wave  statistics  through  regional  ensemble  averaging  and  dense  spatial  coverage  of  both 
wave  and  current  observations. 

To  explore  the  use  of  large  drifter  arrays  to  study  inhomogeneous  wave  conditions 
in  coastal  inlets,  30  WRD  were  deployed  in  the  San  Francisco  Bight.  The  time  series 
recorded  by  free-floating  drifters  in  this  region  with  inhomogeneous  wave  fields  and 
strong  currents  are  generally  nonstationary.  The  analysis  used  here  is  based  on  short  time 
series  to  ensure  statistical  stationarity  so  that  standard  frequency  and  directional  spectral 
analysis  techniques  can  be  used,  and  the  degrees  of  freedom  of  statistical  estimates  are 
increased  through  ensemble  averaging  of  drifter  observations  onto  a  spatial  grid.  The 
resulting  gridded  statistics  provide  a  synoptic  view  of  the  spatial  variability  in  the  wave 
field  and  identify  regions  of  similar  H ^  in  which  the  individual  spectra  estimated  from 
relatively  short  segments  of  the  time  series  were  ensemble  averaged  to  produce  more 
statistically  robust  regional  ensemble  spectra.  These  spectra  were  used  in  conjunction 
with  SWAN  simulations  using  various  combinations  of  physical  processes  and  ray 
computations  (both  with  and  without  currents)  to  identify  the  dominant  physical 
mechanisms  that  result  in  the  spatial  variability  of  the  wave  conditions  in  the  tidal  inlet 
during  strong  ebb  tidal  flow. 

The  spatial  variability  of  the  wave  field  in  the  San  Francisco  Bight  observed  on 
the  April  27,  2012  peak  ebb  tide  was  evaluated  in  the  swell  band,  wind-sea  band  and  the 
high-frequency  tail.  The  variability  in  the  swell  band  was  shown  to  result  primarily  from 
shoaling  and  refractive  effects  as  the  waves  propagate  over  the  San  Francisco  Bar  from 
the  deep  water  of  the  Pacific  Ocean.  East  of  the  bar,  the  combined  refractive  effects  of  the 
shoals  on  either  side  of  the  shipping  channel  results  in  a  divergent  swell  pattern  and  a 
dramatic  reduction  in  swell  wave  energy.  Progressing  farther  east  toward  the  inlet,  the 
dominant  physical  mechanism  influencing  the  wave  field  was  the  ebb  tidal  jet  causing  an 
increase  in  energy  in  the  wind-wave  band  as  these  short  waves  steepened  (mean 
steepness  parameter  increase  from  0.04  to  0.15)  and  refracted  onto  the  ebb  tidal  jet 

(maximum  current  approaching  2  m  s”’ ).  In  the  inlet,  there  is  a  drastic  decrease  in  energy 
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across  the  spectrum  that  can  be  attributed  to  several  factors,  including  depth  induced 
refractive  effects,  the  blocking  and  dissipation  of  short  waves  (frequencies  exceeding 
0.42  Hz)  due  to  the  strong  opposing  currents,  and  a  loss  of  wind  forcing  owing  to 
sheltering  effects  of  the  orography  bordering  the  narrow  inlet. 

These  results  show  that  the  drifter  ensemble  data  resolve  the  two-dimensional 
wave  field  evolution  across  complex  bathymetry  (San  Francisco  Bar)  and  through  an  ebb 
tidal  jet,  providing  a  unique  and  synoptic  analysis  of  the  various  processes  affecting  wave 
dynamics  in  this  complicated  region,  and  illustrating  the  utility  of  this  approach  to  the 
study  of  coastal  inlets. 
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III.  WAVE  EVOLUTION  IN  THE  MOUTH  OF  THE  COLUMBIA 

RIVER 


A.  INTRODUCTION 

Tidal  inlets  and  river  mouths  ean  present  a  significant  hazard  to  navigation  where 
extreme  waves  may  develop  as  waves  interact  with  complicated  shoal  systems  and  strong 
ebb  tidal  jets.  Waves  that  interact  with  currents  and  bathymetry  have  been  studied 
independently.  However,  due  to  the  complications  of  making  measurements  in  such 
dynamic  regions,  little  is  known  of  the  combined  effects  on  the  spatial  wave  evolution. 

As  ocean  waves  transit  regions  with  variable  bathymetry  and  interact  with  strong 
opposing  currents,  they  can  undergo  radical  changes  in  their  physical  properties.  The 
directions  change  due  to  refraction,  wave  heights  increase  as  the  group  speed  slows,  and 
the  waves  can  dramatically  steepen  against  the  current.  In  stationary  conditions,  the 

presence  of  a  current  changes  the  linear  dispersion  relation  to 

1 

ft)  =  U'k  +  (g|k|tanh(|k|/z))^  where  k  is  the  wavenumber  vector,  g  is  gravity  and  the 

absolute  frequency,  ft) ,  is  conserved.  If  the  local  current,  U  ,  and  depth  h  are  assumed 

to  vary  gradually  over  many  wavelengths  (Bretherton  and  Garrett  1968),  it  follows  that  in 

1 

a  framework  that  moves  with  the  medium,  the  relative  frequency,  (7  =  (g|k|tanh(|k|/z))^ , 

experiences  a  Doppler  shift,  cr  =  ft)  -  U  •  k .  The  wavenumber  vector,  k  ,  varies  spatially 
as  well. 

Jonsson  (1990)  provides  a  detailed  overview  of  the  dynamic  and  kinematic  effects 
of  wave-current  interaction.  The  influence  of  currents  on  wave  kinematics  (relative  phase 
velocity  and  wavenumber)  is  highly  dependent  on  the  strength  and  shear  of  the  current 
and  the  angle  at  which  the  waves  approach.  Variations  in  phase  velocity  result  in 
refraction  of  the  waves  that  can  be  described  using  geometric  optics  theory  that  describes 
a  ray  path  normal  to  a  plane  wave,  which  bends  in  response  to  lateral  changes  in  phase 
speed  (Jonsson  1990;  Holthuijsen  and  Tolman  1991).  When  the  waves  oppose  a  strong 
current,  the  reduction  in  group  velocity  and  an  increase  in  relative  frequency  causes  an 
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increase  in  wave  height,  whieh  is  accompanied  by  a  decrease  in  wavelength  that  leads  to 
a  dramatie  steepening  of  the  wave.  These  effeets  are  described  by  the  wave  action 
balance  equation  (eq.  1).  From  the  linear  dispersion  relation,  it  ean  be  shown  that  in 
extreme  eonditions  where  the  waves  direetly  oppose  the  eurrent,  blocking  can  occur 

do 

when  the  relative  group  speed,  c  =  —  equals  the  opposing  eurrent  speed,  U.  The 

ak 

frequeney  at  which  this  occurs  is  known  as  the  blocking  frequency  in  deep  water  given 


by,  fbi 


g 


AnU 


.  Wave  energy  at  frequeneies  greater  than  can  no  longer  propagate. 


Drastic  changes  in  the  wave  conditions  occur  in  blocking  regions  as  the  strong 
nonlinearity  of  waves  breaking  results  in  energy  transfers  to  higher  and  lower  frequeneies 
of  which  only  the  lower  frequency  wave  groups  ean  propagate  against  the  current  (Lai  et 
al.  1989;  Chawla  and  Kirby  2002;  Bascheek  et  al.  2006). 


Wave-current  interaction  in  a  strongly  sheared  eurrent  ean  yield  extreme  wave 
conditions  owing  to  wave  amplifieation  in  eausties  (rays  erossing)  and  the  subsequent 
wave  turning  into  the  eurrent  and  evolving  into  steep  breaking  waves.  These  effeets  were 
observed  as  waves  from  the  southern  oeean  interact  with  the  strong  Agulhus  current 
south  of  the  African  Cape  where  the  Agulhas  current  can  exceed  2  m  s”’  (Irvine  and 
Tilley  1988).  Similar  conditions  develop  in  the  Gulf  Stream  when  waves  oppose  the 
strong  current  (Holthuisen  and  Tolman  1991).  Laboratory  experiments  suggest  that 
enhanced  wave  conditions  can  occur  in  these  situations  partieularly  when  wave-eurrent 
interaction  results  in  refractive  focusing  (Lai  et  al.  1989).  Further,  Chawla  and  Kirby 
(2002)  confirmed  findings  in  Lai  et  al.  (1989)  that  breaking  in  response  to  currents  is  not 
saturated  as  is  found  in  depth-limited  breaking  (Battjes  and  Janssen  1978)  and  that  a 
down  shift  in  the  speetral  energy  oeeurs  as  the  waves  spill  into  the  eurrent.  Additional 
observations  from  field  studies  (Wolf  and  Prandle  1999;  Rusu  et  al.  2011;  Van  der 
Westhuysen  2012;  and  Dodet  et  al.  2013)  outline  the  importance  of  wave-eurrent 
interaetion  on  wave  predietions,  but  consist  of  limited  observations  that  were  made  using 
fixed  instruments. 
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In  tidal  inlets  and  river  mouths  where  complieated  wave  refraetion  patterns  result 
from  the  interaction  with  ebb  tidal  shoal  systems,  observations  are  extremely  limited 
owing  to  the  inherent  difficulties  in  placing  moored  instruments  in  strong  currents,  large 
waves,  and  shipping  lanes  (Thomson  2012;  Elias  et  al.  2012;  and  Pearman  et  al.  2014).  In 
regions  with  a  shoal  system  containing  an  elevated  lens-shaped  bathymetric  feature,  a 
focal  zone,  or  caustic,  develops  in  response  to  rays  converging  as  waves  refract  toward 
shallower  depths  (i.e.  regions  with  slower  group  speeds)  (Vincent  and  Brigg,  1989;  Smit 
and  Janssen,  2013).  It  has  been  theorized  that  nonlinear  instabilities  can  result  in  extreme 
waves  in  these  caustics  (Janssen  and  Herbers  2009).  The  location  of  the  caustic  is  highly 
dependent  on  the  offshore  wave  direction,  which  makes  observing  these  phenomena  with 
traditional  moored  (Eulerian)  sensors  difficult. 

The  mouth  of  the  Columbia  River  (MCR),  known  as  the  ‘Graveyard  of  the  Pacific 
(Oregon  State  Parks),  is  one  of  the  most  dynamic  wave  environments  in  the  world.  It  has 
claimed  approximately  2000  ships  since  1792,  and  serves  as  an  ideal  natural  laboratory 
for  measuring  the  effects  of  a  complicated  shoal  system  and  strong  wave-opposing 
current.  Eai  and  Delisi  (2010)  used  RADAR  observations  to  determine  that  a  large 
degree  of  spatial  variability  exists  in  the  swell  waves  in  the  offshore  approaches  to  the 
MCR.  A  field  study  called  the  Mega-Transect  Experiment  (MGT)  that  was  conducted  by 
the  United  States  Army  Corps  of  Engineers  (US ACE)  in  collaboration  with  the  United 
States  Geological  Service  (USGS)  during  the  summer  of  2005  designed  to  assess  model 
skill  in  the  region  showed  that  wave  current  interactions  are  very  important  as  wave 
heights  nearly  doubled  at  a  station  over  the  bar  (Elias  et  al.  2012).  However,  a  detailed 
spatial  analysis  of  the  physical  mechanisms  that  affect  specific  frequency  band  evolution 
does  not  exist. 

In  this  research,  a  Eagrangian  sensing  strategy  is  employed  in  which  a  large 
number  of  recently  developed  drifters  are  deployed  on  the  peak  of  the  ebb  tidal  current 
and  allowed  to  drift  through  the  most  energetic  wave  regions.  The  spatial  coverage  of  in 
situ  observations  that  results  is  not  practical  using  traditional  Eulerian  (moored)  sensing 
techniques  (Herbers  et  al.  2012,  Thomson  2012,  and  Pearman  et  al.  2014). 
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This  chapter  begins  with  a  deseription  of  the  field  experiment  in  Section  B,  and 
the  analysis  method  is  described  in  Seetion  C.  The  opposing  wave-eurrent  ease  study  is 
presented  in  Seetion  D,  followed  by  a  comparison  of  the  observations  with  SWAN 
predictions  in  Section  E.  Finally,  the  extreme  wave  observations  are  diseussed  in  Seetion 
F,  followed  by  a  summary  in  Seetion  G. 

B,  FIELD  EXPERIMENT 

1.  Field  Site 

The  mouth  of  the  Columbia  River  represents  the  interfaee  between  a  massive  river 
system  and  the  Paeifie  Oeean  (Figure  13).  To  reduee  the  maintenance  requirements  for 
the  navigation  ehannel  and  reduee  the  transit  risk  aeross  the  bar,  the  eonstruetion  of  the 
three-jetty  system  was  eompleted  in  1939  (Figure  13).  A  jet-like  eurrent  develops 
between  the  jetties  as  the  ebb  tide  eouples  with  the  natural  diseharge  of  the  Columbia 
River  that  approaches  8000  m  s  (USGS).  Farge  volumes  of  sediment  are  deposited  as 
the  strong  current  seulpts  the  eomplieated  shoal  system  that  is  known  to  mariners  as  “the 
Bar”.  During  our  deployment,  the  offshore  was  1.75  m  and  the  dominant  offshore 
waves  were  from  the  west-northwest  at  14  seeonds  (Figure  14)  with  seeondary  waves 
from  the  west-northwest  at  10  seeonds  and  a  long  period  (18  seeonds)  south  swell.  In 
addition  to  the  swell,  a  weak  loeally  generated  5-seeond  period  wind-sea  from  the  west- 
northwest  was  present. 

2,  Wave  Resolving  Drifters 

The  three  slightly  different  Wave  Resolving  Drifters  (WRDs)  used  to  measure  the 
eurrents  and  waves  in  this  researeh  (Figure  15)  are  equipped  with  GPS  data  loggers  that 
measure  position  and  horizontal  orbital  veloeities  of  the  swell  waves.  Aeeelerometers 
measure  vertieal  aeeelerations  of  the  higher  frequeney  wind-waves. 

The  first  type  of  WRD,  WRD-A  has  a  1  Hz  GPS  data  logger  that  reeords  u  and  v 
Doppler  veloeity  data  and  positional  data  and  an  aeeelerometer  tethered  beneath  the 
surfaee  drifter  (for  details  see  Meintyre,  2013  and  Pearman  et  al.  2014).  In  the  seeond 
WRD  design,  WRD-B,  an  aerylie  shell  houses  the  same  GPS  data  logger  and  a  more 
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sophisticated  Inertial  Motion  Unit  (IMU)  aoeelerometer  package  (see  Portell,  2013  for 
details).  Both  WRD-A  and  WRD-B  are  also  equipped  with  a  GPS  RF  transmitter  for  real¬ 
time  tracking.  The  final  WRD  design,  WRD-J  is  more  sophisticated  with  a  5  Hz  GPS,  an 
IMU  aeeelerometer  and  an  iridium  transmitter  for  over  the  horizon  traeking,  all  housed  in 
the  aerylic  shell.  The  WRD  are  spherieal  and  have  an  immersion  resonanee  frequeney 
greater  than  1  Hz  that  is  outside  the  range  of  interest  for  this  researeh.  Eaeh  individual 
sensor  has  its  own  time  keeping  method  so  eommon  time  stamps  for  eaeh  drifter  were 
arrived  at  using  the  GPS  time  on  each  WRD  following  the  methods  outlined  in  Herbers  et 
al.  (2012)  and  Pearman  et  al.  (2014). 

3,  Drifter  Deployment 

This  field  experiment  consisted  of  many  deployments;  however,  this  researeh  will 
foeus  on  one  major  deployment  in  whieh  the  offshore  conditions  resulted  in  the  most 
dramatic  wave  evolution.  The  deployment  of  30  WRD  was  eonducted  on  June  08,  2013 
during  the  peak  ebb  tidal  eurrent  from  RW  Point  Sur  as  it  maintained  station  eentered  on 
the  ebb  tidal  jet  in  the  Columbia  River  estuary  (Figure  13),  starting  at  1220  UTC  near  the 
peak  ebb  tidal  eurrent  (predieted  to  oeeur  at  1322  UTC,  NOAA  Tide  Tables). 

The  30  WRD  were  deployed  at  20-seeond  intervals  resulting  in  approximately  60 
m  spread  between  drifters.  As  the  WRD  exited  between  the  jetties  and  transited  the  bar 

with  the  ebb  tidal  jet  exeeeding  veloeities  of  3  m  s”’ ,  they  expanded  laterally  relative  to 
the  mean  drift  direetion  but  remained  elustered  along-traek  (east-west)  following  the 
MCR  alluvial  plume  (see  Figure  16). 

C.  LAGRANGIAN  ANALYSIS 

1.  Drifter  Analysis 

The  drifter  time  series  are  nonstationary  along  the  long  drift  traek  owing  to  the 
Fagrangian  nature  of  the  measurements  as  the  drifters  move  with  the  strong  eurrent 
through  the  spatially  inhomogeneous  wave  field.  Therefore,  two  approaehes  were  taken 
to  the  speetral  analysis.  First,  the  auto-  and  eross-  speetral  estimates  were  eomputed  for 
relatively  short  segments  (204.8-seeonds)  of  the  GPS  and  aeeelerometer  time  series 
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applying  a  Hamming  window.  The  auto  spectra  of  GPS  and  accelerometer  were 
combined  to  form  a  surface  height  spectrum  E{f)  using  the  horizontal  GPS  velocity 
observations  in  the  swell  range  0.03  Hz  -  0.15  Hz  and  vertical  accelerations  in  the  wind- 
sea  range  0.25  Hz  -  1.0  Hz.  In  the  transition  range  of  0.15  Hz  -  0.25  Hz,  a  weighted 
average  of  the  two  sensors  was  used  where  the  GPS  (accelerometer)  weighting  was 
varied  linearly  from  1(0)  at  0.15  Hz  to  0(1)  at  0.25.  Linear  theory  transfer  functions 
(Table  2)  were  applied  to  convert  velocity  and  acceleration  spectra  to  surface  height 

spectra  (see  Pearman  et  ah,  2014  for  details).  From  these  raw  surface  height  spectra  the 

1 

/j.l.0Hz  W 

significant  wave  height,  =  4M  E(f)df\  ,  was  evaluated  (following  Pearman  et 

ah,  2014),  and  these  were  ensemble  averaged  across  drifters  on  a  grid  with  spacing  of  50 
m  (north)  x  100  m  (east)  chosen  based  on  the  lateral  divergence  (north-south)  of  the 
various  drift  trajectories.  To  estimate  accurate  wave  spectra  that  resolve  the  multi-modal 
wave  field,  seven  regions  of  similar  characteristics  were  selected  starting  with  the  Bar  (B) 
and  moving  outside  (inside)  of  the  bar  where  01  (H)  is  closest  and  03  (13)  is  the  farthest 
(Figures  16,  17).  For  each  drifter,  the  entire  time  series  during  passage  of  a  given  region 
was  used  to  estimate  the  surface  height  spectrum,  E(f)  with  FFT  window  lengths  of 

2048  samples  (at  10  Hz)  and  50%  overlap  yielding  a  frequency  resolution  of  0.0049  Hz. 
The  spectrum  from  each  drifter  in  a  given  zone  is  then  ensemble  averaged  and  a  single 
more  statistically  stable  ensemble  spectrum  is  presented  (Figure  18,  discussed  below). 
Additionally,  histograms  of  H  ^  in  each  zone  are  used  to  identify  deviations  from  the 
offshore  distribution  that  can  potentially  indicate  rogue  wave  development  (Figure  19, 
discussed  below). 

2.  Ray-Trajectories 

Rays  based  on  the  offshore  wave  conditions  were  traced  across  the  domain  to 
assist  in  analysis  of  the  wave  field  and,  in  particular,  the  evolution  of  those  waves 
dominant  in  the  spectrum  (see  Figure  20).  Geometric  optics  theory  suggests  that  in  a 
slowly  varying  medium,  a  plane  wave  of  a  particular  frequency  will  follow  the 
wavenumber  vector  that  varies  slowly  in  response  to  refraction  by  depth  and 
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current  variations.  The  wave  number  magnitude  k  is  ealeulated  from  the  dispersion 
relation  based  on  the  relative  frequeney  G  and  loeal  depth.  The  rays  are  traeed  with  a 
simple  iterative  seheme  to  adjust  (7  for  the  Doppler  shift  k  •  U  that  depends  on  the  angle 
of  the  loeal  plane  wave  to  the  loeal  eurrent  from  the  modeled  flow  field  deseribed  in 
Seetion  E.  In  eaeh  iteration  the  ray  equation: 


ds 


- - sin(0)  — +  cos(0)  — 

rcQ  ^  UX  uy 


(3) 


is  solved  where  — are  ealeulated  on  a  50  m  x  50  m  grid.  The  starting  points  for  the 
ax  ay 


rays  are  seleeted  at  an  offshore  loeation  sueh  that  the  bathymetry  and  eurrents  do  not 
influenee  the  waves.  The  ray  spaeing  is  eonstant  along  the  outer  edge  of  the  grid  (50  m) 
and  eaeh  ray  is  traeed  using  10  m  inerements  along  the  ray  trajeetory,  s  ,  and  analyzing 


using  aneillary  relations 


dx 

ds 


=  COS0, 


—  =  sin0 . 
ds 


The  predieted  bend  in  response  to  the  phase  speed  variations  along  the  wave  erest 
indueed  by  the  variable  bathymetry  and  sheared  eurrents.  The  eonvergenee  (divergenee) 
of  the  rays  that  results  indieates  areas  of  inereasing  (deereasing)  energy.  When  the 
refraetion  is  strong  enough,  a  eaustie  (rays  erossing)  ean  result,  and  when  eurrent  shear  is 
strong  the  rays  ean  beeome  trapped  and  refraet  baek  and  forth  aeross  the  eurrent  in  a 
snake-like  pattern  (Mei  1983). 


D,  CASE  STUDY:  WAVES  OPPOSING  CURRENT 

The  offshore  wave  eonditions  are  somewhat  eomplex;  however,  the  following 
analysis  will  foeus  on  the  dominant  west-northwest  swell  (14-seeond)  and  the  5- 
seeond  wind-sea  that  are  affeeted  in  different  ways  by  the  shoal  and  tidal  eurrents  that 
exeeed  3  m  s”'  (Figure  16).  This  approaeh  eapitalizes  on  the  dense  spatial  eoverage 
of  H ^  to  obtain  a  synoptie  view  of  the  wave  height  variability  aeross  the  mouth  of  the 
Columbia  River  (Figure  17).  Additionally  wave  speetral  estimates  in  the  seven 
seleeted  zones  (Figure  18)  are  presented  to  provide  a  more  detailed  analysis  of  the 
transformation  of  the  swell  and  sea  eomponents.  The  wave  height  and  speetral 
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estimates  are  eombined  with  ray  eomputations  to  identify  the  bottom  and  eurrent 
features  that  are  primarily  responsible  for  the  wave  variability. 

1.  Dominant  Wave  Evolution 

The  offshore  H ^  are  relatively  homogenous  with  a  mean  H ^  near  2.5  m  and  a 
standard  deviation  in  the  distribution  of  about  0.9  m  (Figure  19).  The  peak  of  the  speetrum 
(0.065  Hz)  in  Zone  03  is  dominated  by  a  14-seeond  swell  where  the  direetional  spread  is  at  a 
minimum  and  the  direetion  is  from  the  west-northwest  (Figure  18).  The  straight  and  evenly 
spaeed  rays  indieate  little  bathymetrie  or  eurrent  influenee  (Figure  20)  until  the  lens  shaped 
shoal  in  the  NE  eomer  of  Zone  02  where  the  rays  refraet  and  eonverge  in  a  well-defined 
eaustie.  The  resulting  eross-wave  patterns  are  evident  in  the  nearby  doubling  of  the 
direetional  spread  (Figure  18).  The  spread  in  H ^  estimates  remains  near  eonstant  although 
the  mean  inereases  to  2.8  m.  The  region  of  inereased  H ^  in  the  eenter  of  the  zone  (Figure  17) 
elosely  eoineides  with  the  eore  of  the  ebb  tidal  jet  (Figure  16)  suggesting  some  shoaling  (i.e. 
deereasing  group  speed)  in  response  to  the  opposing  eurrent  (>  2.5  m  s”' )  may  also 
eontribute  to  the  enhaneed  wave  heights.  This  wave-eurrent  interaetion  is  also  evident  in  the 
Doppler  shift  of  the  speetral  peak  to  0.083  Hz. 

In  Zone  01  the  elevated  peak  in  the  speetrum  is  further  Doppler  shifted  (Figure 
18)  as  the  opposing  eurrent  reaehes  3  m  s”'  (Figure  16),  and  a  foeused  line  of  elevated 
oriented  from  the  west-northwest  to  the  east-southeast  develops  (Figure  17)  that 
elosely  eorresponds  to  the  eaustie  that  forms  as  the  rays  eonverge  down  wave  of  the  lens 
shaped  shoal  loeated  in  Zone  02  (Figure  20).  The  enhaneed  H ^  of  3.7  m  are  almost 
double  the  ^  2  m  H ^  values  observed  to  the  north  and  south  of  the  foeusing  region. 
Similar  enhaneements  were  identified  in  laboratory  experiments  (Vineent  and  Briggs, 
1989)  where  narrow  band  (in  a  direetional  sense)  waves  were  found  to  nearly  double  in 
height  in  the  eausties. 

In  Zone  B,  on  the  bar,  the  wave  heights  are  further  enhaneed  (H ^  >  5  m),  with  the 
largest  waves  in  the  northern  portion  of  this  zone  in  eontrast  to  the  -'3  m  observed  in 
southern  portion.  (See  also  the  signifieantly  greater  spread  in  the  distribution.  Figure  19) 
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The  sharp  gradient  of  wave  heights  over  very  short  distanees  (Figure  17)  suggests  that 
these  enhaneed  waves  are  the  result  of  loeal  physieal  proeesses.  The  inerease  in  energy 
appears  to  be  concentrated  at  the  peak  of  the  spectrum  (Figure  18).  However,  the  rays 
from  the  dominant  280  degrees  offshore  direction  pass  to  the  south  of  the  region  of 
maximum  wave  height.  The  refraction  effects  are  sensitive  to  the  incident  wave  direction, 
and  therefore,  it  is  important  to  consider  the  directional  spread  about  the  mean  wave 
direction.  The  directional  spread  is  15  degrees  near  the  spectral  peak  (Figure  18,  bottom), 
so  waves  approach  the  shoal  from  a  range  of  angles  around  the  dominant  offshore 
direction.  Rays  at  the  same  offshore  frequency  from  270  degrees  (Figure  21)  refract 
around  the  south  side  of  the  shoal  and  are  indeed  concentrated  in  this  region  of  enhanced 
H ^  well  down-wave  from  the  caustic  providing  a  possible  explanation  for  this 

pronounced  enhancement. 

Behind  the  bar  in  Zone  11  the  enhanced  H ^  are  no  longer  evident,  which  is 
consistent  with  the  reduction  in  the  peak  of  the  spectrum  (Figure  18).  The  persistent 
strong  current  and  lack  of  a  noticeable  Doppler  shift  suggests  that  the  notable  decrease  in 
H ^  on  the  central  portion  of  the  zone  can  be  partially  attributed  to  the  increased  water 
depth  (unshoaling),  but  more  importantly,  the  ray  trajectories  in  this  deeper  area 
surrounded  by  shoals  indicate  an  energy  divergence  similar  to  observations  in  other  inlet 
studies  (Elias  et  al.  2012;  Dodet  et  al.  2012;  Pearman  et  al.  2014). 

Zone  12  is  well  within  the  inlet,  however,  H ^  remains  elevated  near  3  m  (Figure 
17),  and  the  peak  shifts  further  to  a  higher  frequency  (Figure  18).  The  overall  pattern 
shows  ray  divergence  toward  the  sides  of  the  inlet,  and  a  concentration  of  the  rays  in  this 
region  (Figures  20  and  21)  turns  into  the  current  suggesting  waves  in  this  band  continue 
to  shoal  into  the  strong  current  (>3  ms”’).  The  change  in  the  mean  wave  direction  from 
southwest  to  west  and  low  directional  spread  (Figure  18)  are  consistent  with  this  ray 
pattern. 

By  the  time  the  dominant  rays  reach  the  eastern  portion  of  Zone  13  the  divergence 
of  the  wave  energy  results  in  a  dramatic  decrease  in  spectral  levels  at  the  dominant  swell 
peak  (Figure  18).  The  decrease  in  H ^  to  about  1  m  (Figure  17)  is  similar  to  what  was 


27 


observed  in  SF  Bay  (Pearman  et  al.,  2014).  However,  it  occurs  much  farther  into  the  inlet 
in  MCR  where  the  currents  are  significantly  stronger  (about  a  factor  of  2)  indicating  that 
wave  focusing  on  the  current,  which  counteracts  energy  divergence  caused  by  depth- 
induced  refraction,  is  more  important  in  the  MCR  than  in  SF  Bay. 

2,  Wind-wave  Band 

Interpretation  of  the  evolution  of  the  wind-wave  band  is  complicated  by  the 
directional  spread,  which  is  much  larger  than  at  the  dominant  peak  (Figure  18). 
Laboratory  experiments  (Vincent  and  Briggs,  1989)  show  that  for  broadband  (directional) 
waves  the  increase  in  wave  heights  on  opposing  currents  is  much  less  dramatic  than 
narrow-band  waves.  Therefore,  when  calculating  the  blocking  frequency  (the  frequency 
at  which  waves  will  no  longer  propagate  into  an  opposing  current)  using  the  linear 
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dispersion  relation,  = - ,  where  a  is  the  difference  in  angle  between  -  U  and 

k ,  it  is  important  to  note  that  only  the  waves  for  which  a  is  sufficiently  small  will 
actually  be  blocked  (e.g.  is  twice  as  large  for  an  ct  =  60°  compared  with  the  pure 
opposing  cc  =  0  case).  Thus,  when  a  directionally  broad  wave  field  encounters  a  current 
jet,  the  small  (X  components  are  amplified  as  they  refract  and  shoal  into  the  current, 
whereas  the  large  (X  components  merely  cross  the  jet.  Bathymetry  effects  are  negligible 
as  is  clear  in  the  comparison  of  the  rays  with  currents  (Figure  22)  and  without  currents 
(Figure  23).  The  following  analysis  progresses  from  the  offshore  Zone  03  to  inside  of 
the  bar  to  in  Zone  13. 

The  relatively  broad  spectral  peak  near  0.2  Hz  in  Zone  03  is  amplified 
progressively  until  it  reaches  a  maximum  in  Zone  B  where  the  mean  velocity  is  3  m  s“' 
and  f|^^  =  021  Hz  for  a  =  0 .  The  considerable  convergence  of  the  rays,  particularly  in 
the  eastern  portion  of  the  zone,  indicates  that  the  waves  are  becoming  highly  focused  on 
the  current  and  thus  refraction  (in  addition  to  shoaling)  may  contribute  to  the  increase  in 
energy.  The  spectrum  decreases  slightly  in  Zone  11  before  increasing  in  Zone  12  where 
the  rays  are  now  directly  opposing  the  current,  and  the  directional  spread  has  reduced 
from  30  degrees  to  20  degrees.  The  dramatic  reduction  in  energy  in  Zone  13  is  indicative 
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of  blocking  of  those  waves  angled  into  the  eurrent  in  Zone  12  where  the  mean  veloeity  is 
3.3  m  s“'  and  the  =  0.23  Hz  for  a  =  0  .  Additionally,  the  ray  trajeetories  indieate  that 
waves  from  the  mean  offshore  direetion  beeome  trapped  (refraeting  baek  and  forth  aeross 
the  jet  as  deseribed  in  Mel  1983)  on  the  eurrent  which  by-passed  Zone  13  almost  entirely 
(Figure  22). 

E.  SWAN  COMPARISON 

The  wave  aetion  balanee  equation  at  the  core  of  phase  averaged  models  like 
SWAN  has  been  shown  to  prediet,  with  reasonable  skill,  the  evolving  wave  statistics  over 
slowly  varying  bathymetry  and  in  the  presenee  of  currents  (Booij  et  al.  1999  and  Ris  et  al. 
1999).  However,  phase-averaged  models  tend  to  generally  over-prediet  in  strong 
opposing  eurrents  (Van  der  Westhuysen  2007  and  2012)  and  to  miss  the  peaks  in  regions 
of  highly  foeused  wave  energy  (Cavaleri  2009;  Janssen  and  Herbers  2009;  Smit  and 
Janssen  2013).  The  measurements  used  in  prior  field  evaluations  of  SWAN  were  mostly 
Eulerian  in  nature,  relatively  sparse  and  do  not  cover  wave  evolution  in  extreme  currents 
and  through  caustic  regions. 

In  this  researeh  SWAN  was  one-way  eoupled  with  veloeity  output  from  the 
eireulation  model,  SELFE  (semi-implieit  Eulerian-Eagrangian  finite-element  model). 
SEEFE  is  a  4-dimensional  baroelinio  eireulation  model  that  has  been  designed  for  use  in 
the  MCR  region  (Zhang  and  Baptista  2008).  SEEFE  solves  the  3D  nonlinear  shallow- 
water  equations  for  free  surface  elevation  and  3D  veloeity-salinity-temperature  with 
hydrostatie  and  Boussinesq  approximations  and  transport  equations  for  salt  and  heat.  The 
external  foreing  eonsists  of  4-semidurnal  tidal  eonstituents  along  the  offshore  boundary 
and  the  diseharge  of  the  Columbia  River  at  a  loeation  suffieiently  upstream  to  be  beyond 
the  maximum  salt-water  intrusion,  using  measurements  gathered  as  part  of  from  Center 
of  Coastal  Margin  Observation  and  Predietion  (CMOP)  program.  (SEEFE  model  output 
is  downloadable  from  the  CMOP  website:  http://www.steemop.org) 

The  SWAN  model  was  implemented  on  a  variable  resolution  grid  that  uses  -'100 
m  grid  spaeing  offshore  and  10  m  in  regions  where  bathymetry  (from  the  NOAA 
bathymetry  database)  showed  rapid  horizontal  variability.  The  wind  and  wave  foreing 
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were  taken  from  observations  made  at  NOAA  buoy  46089  loeated  160  km  offshore  (for 
details  on  souree  term  seleetion  and  set  up  see  Pearman  et  al.  2014).  Multiple  SWAN 
simulations  in  whieh  souree  terms  were  individually  disabled  were  eondueted  to  evaluate 
the  dominant  physieal  meehanisms  influeneing  the  SWAN  predietions  (Table  1).  The 
from  the  SWAN  simulations  (Figure  24)  evaluated  here  are  eolloeated  in  spaee  and  time 
with  WRD  measurements  (Figure  17),  and  the  SWAN  speetra  (Figure  25)  are  seleeted  at 
a  grid  point  elose  to  the  eenter  of  the  eaeh  zone  at  a  time  eonsistent  with  the  drifter 
observations  from  that  zone. 

The  general  wave  pattern  predieted  by  the  SWAN  simulation  (Figure  24) 
eompares  rather  well  with  the  WRD  estimates  (Figure  17).  However,  the  SWAN 
predietions  tend  to  be  biased  high,  partieularly  in  Zones  12  and  13  where  the  eurrents  are 
strongest  (Figure  26).  Of  partieular  interest  is  the  elear  jump  in  the  peak  of  the  speetrum 
moving  from  Zone  02  to  Zone  01  in  response  to  the  eaustie  that  forms  down- wave  from 
the  lens-shaped  shoal  (Figure  13).  The  loeation  of  the  elevated  H^in  the  SWAN 
simulation  (Figure  12)  eoineides  with  the  loeation  where  the  WRD  observe  amplifieation, 
but  show  about  a  +0.5  m  high  bias.  The  elevated  energy  level  at  the  peak  of  the  speetrum 
when  eompared  to  the  WRD  in  this  zone  may  be  due  to  the  ensemble  averaging  of  the 
WRD  speetrum  ealeulated  within  this  zone  versus  the  single  point  speetrum  of  SWAN 
that  falls  within  the  eaustie  in  the  ray  trajeetories  (Figure  21). 

The  region  of  extreme  that  was  observed  in  the  WRD  measurements  in  Zone 

B  (Figure  17)  is  evident  in  the  SWAN  predietions  slightly  further  east  (on  the  edge  of 
Zone  B  and  Zone  11).  However,  the  enhaneed  are  under  predieted  by  as  mueh  as  3.5 
m  or  50%,  and  the  high  degree  of  spatial  variability  in  the  WRD  measurements  when 
eompared  to  the  smooth  SWAN  predietions  demonstrate  that  SWAN  does  not  effeetively 
eapture  this  sharp  region  of  very  enhaneed  waves  (Figure  26).  While  the  peak  of  the 
speetrum  remains  elevated  in  this  zone  (Figure  25)  it  is  fairly  elose  to  the  previous  zone 
and  does  not  show  the  dramatie  inerease  in  the  dominant  speetral  peak  that  is  evident  in 
the  WRD  speetrum  (Figure  18).  Rather  the  limited  inerease  in  results  from  elevated 

energy  levels  at  higher  frequeneies  in  response  to  the  strong  eurrent  (-'3  m  s”’ ),  whieh  is 
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confirmed  by  the  mueh  lower  H ^  from  the  SWAN  simulation  in  whieh  the  eurrents  are 
disabled  (Figures  26  and  27). 

The  next  region  of  elevated  wave  heights  falls  in  Zone  12  where  SWAN  prediets 
H ^  in  excess  of  4m  eompared  to  the  observed  wave  heights  of  less  than  3  m  (Figure  24). 
The  SWAN  simulation  with  the  eurrents  disabled  (Figure  27)  indieates  that  the  currents 
are  responsible  for  the  high  bias,  whieh  is  eonfirmed  in  the  1-1  eomparison  (Figure  26). 
The  speetrum  in  this  zone  (Figure  25)  shows  that  the  elevated  energy  levels  at 
frequeneies  higher  than  the  speetral  peak  and  provides  a  elear  example  of  the  model 
under-dissipating  wave  energy  on  strong  opposing  eurrents. 

F,  EXTREME  WAVES 

Whereas  mueh  of  the  spatial  wave  variability  observed  with  the  WRD  array  is  at 
least  qualitatively  eonsistent  with  expeeted  bottom-  and  eurrent-  induced  shoaling  and 
refraetion  effeets,  the  extreme  waves  observed  in  Zone  B  warrant  more  discussion.  The 
dramatie  changes  in  H ^  that  oeeur  on  relatively  small  spatial  seales  over  the  MCR  Bar 
(zone  B)  are  highly  loealized  in  the  eentral  and  northern  portion  of  the  zone  (Figure  17). 
Multiple  drifters  experienced  enhanced  wave  development  in  a  very  foeused  region  that 
is  well  down-wave  from  the  refraction  caustic  (Figure  28).  The  waves  nearly  triple  in  H ^ 
over  a  distance  of  less  than  a  few  wavelengths.  Comparison  of  the  speetra  that  have  been 
ensemble  averaged  over  enhaneed  >=  4.5  m)  and  unenhaneed  (77^ <4.5  m) 
individual  drifter  time  series  shows  that  the  most  signifieant  amplifieation  in  the  speetrum 
oeeurs  at  the  dominant  peak  (Figure  30). 

The  waves  become  foeused  on  the  eombination  of  the  ridge  shaped  shoal  and  the 
opposing  eurrent  and  the  direetional  spread  is  7  degrees  narrower  at  the  speetral  peak  in 
the  enhaneed  speetrum  than  it  is  in  the  unenhaneed  speetrum  (Figure  30).  A  elear 
negative  eorrelation  between  77^  and  direetional  spread  (Figure  31)  supports  the 
hypothesis  that  refraetive  foeusing  eauses  the  anomalously  larger  waves. 

To  further  investigate  the  behavior  of  these  waves,  it  is  helpful  to  do  so  in  the 
time  domain.  However,  as  outlined  earlier,  the  aeeelerometer  (direet  measure  of  vertieal 
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motions)  is  inaccurate  in  the  swell  band  where  double  integration  of  the  aeeelerations  to 
arrive  at  displaeements  results  in  an  0(4)  amplifieation  of  the  noise  in  the  swell  band  (see 
Appendix  A).  However,  in  a  narrow  band  linear  wave  field  the  horizontal  veloeity 
eomponent  in  the  dominant  wave  direetion  is  proportional  to  the  vertieal  wave  height 
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profile  with  a  transfer  function  that  equals  —  (Herbers  et  al.  2012).  Therefore,  the 
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horizontal  orbital  motions  in  the  time  domain  provide  insight  into  the  wave  evolution 
(Figure  31).  The  drifter  shown  here  that  reeords  the  largest  waves  is  moving  from  right  to 
left,  and  thus  the  evolution  of  waves  traveling  from  left  to  right  is  shown  in  reverse. 
Between  times  12:40  and  12:34  the  orbital  veloeities  are  unremarkable  with  the 
appearanee  of  normally  shaped  wave  groups  and  variations  ±1  m  s~' .  As  the  drifter 
passes  through  the  amplified  wave  group  (near  time  12:34-12:32),  the  orbital  veloeities 
beeome  so  extreme  as  to  oeeasionally  exeeed  the  mean  drift  veloeity  of  3  m  s~'  and  a 
eomplete  reversal  of  flow  results.  These  enhanced  orbital  veloeities  are  more  than  twiee 
those  measured  before  and  after  the  extreme  wave  event.  The  highly  loealized  nature  of 
these  enhaneed  veloeities  suggests  rogue  wave  development. 

The  rays  in  this  region  are  within  a  few  degrees  of  the  primary  eurrent  direetion 
(Figure  28),  eonditions  that  were  found  to  be  ideal  in  deep  water  investigations  of  rogue 
wave  development  diseussed  in  the  introduction.  While  the  lens  shaped  shoal  is  eritical 
in  setting  up  the  conditions  for  rogue  wave  development,  the  actual  extreme  waves  oecur 
well  down  wave  from  the  dramatic  caustic  as  the  waves  beeome  focused  in  response  to 
the  ridge  shaped  shoal  and  the  strong  eurrent.  This  indieates  that  the  development  is  in 
response  to  nonlinear  instabilities  that  steepen  and  enhanee  the  waves  in  this  region, 
similar  to  what  was  modeled  in  the  idealized  eonditions  by  Janssen  and  Herbers  (2009). 

G.  SUMMARY 

The  eomplieated  bathymetry  that  is  seulpted  by  the  strong  outflow  from  the 
Columbia  River  ereates  a  spatially  variable  wave  field  that  is  diffieult  to  measure  using 
traditional  moored  instruments.  The  lens-shaped  shoal  at  the  leading  edge  of  the  bar 
system  is  the  dominant  forcing  feature  in  the  region  that  focuses  swell  waves  through  a 
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caustic  as  they  refract  toward  shallow  water.  The  waves  then  encounter  a  eomplicated  bar 
system  that  has  a  ridge  shaped  shoal  extending  out  the  lens  and  a  shallow  bar  that  crosses 
the  ehannel.  The  effeets  of  this  irregular  bathymetry  and  the  combined  tidal  and  river 
flow  that  ean  exeeed  4  m  s~'  result  in  loealized  extreme  wave  eonditions. 

To  measure  this  wave  evolution  30  Wave  Resolving  Drifters  (WRD)  were 
deployed  on  the  peak  ebb  tide  from  inside  the  river  mouth  and  allowed  to  drift  through 
the  opposing  wave  field  with  the  eurrent  and  reeovered  at  the  seaward  edge  of  the  river 
plume  several  kilometers  offshore.  Beeause  of  the  exeeedingly  strong  current  and  general 
spatial  inhomogeneity  of  the  wave  field  through  whieh  the  WRD  passed,  relatively  short 
segments  were  used  to  estimate  significant  wave  height,  H ^ ,  whieh  were  then  ensemble 
averaged  on  a  spatial  grid  in  order  to  obtain  a  robust  synoptic  estimate  of  the  wave  field 
variability.  Zones  of  similar  wave  eharaeteristies  were  identified  and  speetra  were 
caleulated  from  the  time  series  of  the  drifter  measurements  within  the  zone  and  ensemble 
averaged  to  resolve  the  speetral  properties  of  the  eomplex  multi-modal  wave  field. 

The  offshore  wave  conditions  on  8  June  2013  eonsisted  of  a  dominant  14-second 
swell  and  a  weak  5 -second  wind-sea  both  arriving  from  the  west-northwest.  The  wind-sea 
was  highly  influenced  by  the  strong  ebb  tidal  jet  exiting  the  channel.  It  was  shown  using 
rays  that  the  waves  appear  to  beeome  trapped  on  the  strong  eurrent.  SWAN-SELFE 
coupled  model  simulations  were  eompared  with  these  observations  and  it  was  found  that 
SWAN  simulations  tend  to  under-dissipate  the  high  frequeney  waves  as  they  shoal  and 
break  onto  the  strong  eurrent  and  generally  earry  too  mueh  energy  in  the  higher 
frequencies  well  into  river. 

The  dominant  west-northwest  14-seeond  period  waves  refraet  on  the  lens  shaped 
shoal  at  the  outer  approaches  to  the  MCR  ebb  tidal  shoal  region  and  form  a  eaustie  on  the 
eastern  side.  The  wave  heights  are  amplified  in  the  eaustie  region  largely  due  to  the 
refraetive  effeets.  However,  the  rays  that  refraet  around  the  southern  part  of  the  lens 
remain  tightly  spaeed  as  they  exit  the  eaustie  and  oppose  the  current  while  moving  along 
a  ridge  shaped  shoal  that  extends  seaward  from  the  north  jetty.  Extreme  waves  were 
observed  along  these  rays,  and  although  there  is  enhaneement  in  the  high-frequeney 
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portion  of  the  spectrum,  the  dominant  increase  in  energy  is  within  the  swell  peak.  The 
observed  WRD  orbital  velocities  in  the  region  show  an  abrupt  doubling  in  intensity 
indicating  rogue  wave  development. 

A  coupled  SWAN-SELFE  simulation  shows  a  slight  wave  enhancement  in  the 
region  where  large  waves  were  observed.  However,  the  increased  energy  levels  predicted 
by  SWAN  are  in  the  higher-frequency  portion  of  the  spectrum  (away  from  the  swell 
peak)  indicating  the  model  does  not  adequately  dissipate  waves  that  are  blocked  by  the 
opposing  current.  The  lack  of  enhancement  at  the  spectral  peak  in  SWAN  supports  the 
hypothesis  that  the  rogue  wave  development  is  largely  due  to  nonlinear  instabilities  that 
develop  down  wave  from  the  caustic  rather  than  from  purely  refractive  effects. 
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IV.  SUMMARY  AND  FUTURE  WORK 


A.  SUMMARY 

The  two  chapters  contained  herein  outline  the  development  and  application  of 
new  drifters  that  were  engineered  for  use  in  tidal  inlets  and  estuaries.  The  wave  resolving 
drifters  (WRD)  were  found  to  produce  results  comparable  with  more  sophisticated  and 
expensive  ‘state  of  the  art’  three  axis  wave  buoys.  The  deployment  and  of  the  WRD  in 
large  numbers  coupled  with  a  methodology  based  on  ensemble  statistics  of  the  non¬ 
stationary  time  series  facilitated  evaluation  of  wave  evolution  at  a  spatial  density  that  was 
not  achievable  using  traditional  (Eulerian)  in  situ  measurement  strategies  alone. 

The  addition  of  an  accelerometer  package  removes  the  180-degree  direction 
ambiguity  that  was  evident  in  earlier  GPS-only  drifter  designs.  The  accelerometer  also 
provides  direct  vertical  measurements  of  higher  frequency  wind  waves  in  the  spectral  tail 
that  is  not  possible  using  GPS  that  are  less  responsive  at  high  frequencies.  Ensemble 
averaging  wave  spectra  over  many  drifters  allowed  us  to  resolve  the  complex 
inhomogeneous  and  multi-modal  wave  fields.  These  spectra  were  used  in  conjunction 
with  ray  computations  to  provide  insight  into  the  frequency-specific  evolution  of  a  wave 
field. 

The  physical  mechanisms  responsible  for  the  variable  wave  field  in  the  San 
Erancisco  Bight  were  analyzed  with  coupled  DelftSD  EEOW  and  SWAN  simulations.  It 
was  shown  that  refraction  over  the  bifurcated  shoal  system  of  the  Eour  Eathom  Bank  and 
the  South  Bar  is  the  dominant  mechanism  that  controls  the  spatial  variability  of  the  swell 
waves.  In  contrast,  the  strong  ebb  tidal  jet  controls  the  wind  wave  band.  The  shape  of  the 
inlet  results  in  divergent  rays  at  the  both  the  swell  and  wind  wave  frequencies,  and  thus 
wave  energy  entering  the  channel  is  limited. 

The  conditions  in  the  Mouth  of  the  Columbia  River  (MCR)  were  challenging. 
Both  wave-bottom  and  wave-current  interaction  processes  contributed  to  the  wave 
evolution  similar  to  what  was  observed  in  the  San  Erancisco  Bight.  However,  unlike  the 
San  Erancisco  Bight,  in  the  MCR  the  wind  waves  get  trapped  on  the  current  and  penetrate 
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well  into  the  river  mouth.  Additionally,  the  longer  period  west  swell  is  focused  sharply 
on  the  lens-shaped  shoal  outside  of  the  jetties  forming  a  well-defined  caustic.  The  largest 
waves  were  observed  well  down-wave  from  the  caustic  indicating  likely  rogue  wave 
development  in  response  to  nonlinear  instabilities  of  waves  amplified  by  wave-bottom 
and  wave  current  interaction.  This  experiment  confirmed  that  the  phase-averaged  models 
used  operationally  by  the  Navy  tend  to  under-dissipate  wave  energy  in  strongly  opposing 
currents,  but  it  also  highlighted  that  these  models  do  not  resolve  and,  in  fact,  drastically 
under-predict  extreme  wave  development  in  caustic  regions. 

B,  FUTURE  WORK 

The  analysis  of  these  experiments  is  ongoing.  The  additional  deployments  that 
were  conducted  in  the  San  Francisco  Bay  provide  an  opportunity  to  investigate  how 
waves  develop  in  these  regions  under  various  offshore  conditions.  On  February  15,  2012, 
the  drifters  followed  the  complicated  current  pattern  in  which  some  of  the  drifters  went 
offshore  over  the  bar,  but  one  WRD  drifted  north  over  the  Four  Fathom  Bank,  was 
trapped  in  an  eddy,  crossed  the  bank  again  and  eventually  entered  and  exited  the  surf 
zone.  Another  drifter  exited  the  channel  and  did  not  cross  the  South  Bar,  before  the 
current  reversed  and  eventually  drifted  along  the  coast  of  Ocean  Beach.  Several  large 
deployments  were  conducted  in  October  of  2012  to  observe  wave  evolution  across  the 
bar  in  more  detail. 

Two  additional  major  deployments  were  conducted  at  the  Mouth  of  the  Columbia 
River.  The  first  was  on  May  27,  2013,  when  there  was  a  south  swell  and  the  winds  were 
offshore.  The  preliminary  analysis  of  the  30  WRD  time  series  and  evolution  of  the 
spectrum  indicate  fetch  limited  wave  growth  on  a  strong  following  current.  In  the  second 
major  deployment  on  May  29,  2013,  a  total  of  47  WRD  were  deployed  and  allowed  to 
drift  approximately  18  km  north  following  the  ebb  tidal  plume.  On  both  of  these  days 
extreme  wave  development  was  not  observed,  which  indicates  that  the  offshore  wave 
conditions  may  be  an  important  factor  in  establishing  the  focusing  that  leads  to  rogue 
wave  development. 
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C.  NAVY  RELEVANCE 

Recent  advances  in  GPS  and  Inertial  Measurement  Units  (IMU)  technology  have 
significantly  reduced  the  total  cost  for  the  components  used  in  these  wave-resolving 
drifters  to  around  $1000.  The  processing  that  was  done  here  could  be  completed  onboard, 
thus  allowing  for  spectra  and  bulk  wave  statistics  (a  relatively  small  data  set)  to  be 
transmitted  in  near  real-time.  Further,  the  error  in  the  IMU  measurements  is  significantly 
reduced  when  the  3 -axis  magnetometer,  3 -Axis  gyro  and  3 -axis  accelerometer 
combination  is  used  to  transform  the  body  frame  accelerations  into  linear  accelerations  in 
the  earth  reference  frame  (see  Portell,  2013).  It  follows  that  as  the  technology  (software 
and  hardware)  improves,  the  double  integration  from  accelerations  to  displacements  in 
the  time  domain  will  facilitate  direct  investigation  of  individual  wave  parameters  (slope, 
height,  and  length)  that  are  likely  important  in  understanding  rogue  wave  development. 

Drifters  like  these  can  also  be  used  to  fill  in  the  gaps  in  regions  where  the  existing 
observational  network  is  limited,  or  more  detailed  spatial  resolution  is  necessary  for  the 
validation  and  improvement  of  numerical  model  predictions.  While  a  detailed  model 
validation  is  not  the  focus  of  this  research,  it  clearly  shows  that  these  drifters  are  useful  in 
identifying  the  weaknesses  in  the  phase-averaged  models  that  are  currently  used  by  the 
U.S.  Navy  to  make  predictions  in  shallow  water  regions.  Current  literature  that  identifies 
the  tendency  for  SWAN  to  over-predict  the  wave  heights  in  regions  with  waves  opposing 
strong  currents  is  limited  by  the  Eulerian  observations  that  were  used.  However,  the 
comparison  in  this  research,  made  possible  by  these  spatially  dense  Lagrangian  drifter 
observations,  shows  that  SWAN  also  severely  under  predicts  the  peaks  of  extreme  wave 
amplification  in  these  conditions.  An  assumption  that  waves  are  over  predicted  can  result 
in  elevated  risk  to  the  mission  and  an  increased  risk  to  forces  resulting  from  the  very  high 
probability  that  operators  will  experience  wave  conditions  that  well  exceed  their 
operational  limits.  The  improved  understanding  resulting  from  this  research  can  help 
eliminate  this  worst-case  operational  risk  management  scenario. 
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APPENDIX  A.  VERTICAL  DISPLACEMENT  SPECTRA 
AND  BULK  STATISTICS 

We  evaluated  the  new  WRD  performanee  through  eomparison  with 
measurements  from  a  Datawell  Wave  Rider  DWR-G4  direetional  buoy,  which  derives 
three-dimensional  velocity  from  the  Doppler  shift  estimates  of  the  raw  GPS  LI  carrier 
signal  (similar  to  the  GT-31  GPS).  In  the  internal  processing  these  velocities  are 
integrated  to  produce  three-dimensional  displacement  time  series.  The  Datawell  GPS- 
based  measurements  have  been  carefully  validated  against  earlier  generation  Datawell 
buoys  with  accelerometers,  compass  and  tilt  sensors  (de  Vries  et  al.,  2003,  Herbers  et  al., 
2012)  and  thus  can  be  considered  a  reliable  “ground  truth”  for  the  WRD  measurements. 

We  selected  a  deep-water  site  (1700  m)  well  offshore  (40  km)  to  reduce  effects 
from  coastal  sheltering,  shoaling  and  refraction  on  the  wave  statistics.  A  well-developed 
wind  sea  in  a  nearly  constant  northwest  wind  of  '-'15  m  s~'  persisted  throughout  the 
period  over  which  the  measurements  were  made.  Three  WRDs  were  deployed  alongside 
the  Datawell.  The  largest  separation  distance  between  the  Datawell  and  the  WRDs  during 
the  experiment  was  2.1  km,  which  is  small  compared  with  the  typical  0(100  km)  scale  of 
variability  in  open  ocean  wave  conditions.  Therefore,  the  statistics  of  the  wave  field 
measured  by  the  WRD  and  the  Datawell  buoy  can  be  assumed  to  be  the  same. 

A.  ACCELEROMETER  PERFORMANCE 

The  X6-2  sensor  was  tethered  beneath  the  drifter  to  decouple  the  sensor  from  the 
pitch  and  roll  motions  of  the  surface  buoy  in  steep  waves.  However,  while  the  pitch  and 
roll  effects  are  effectively  dampened,  occasional  steep,  near-breaking  waves  may  cause 
large  tilt  angles  of  the  ballast  chain  and  sensor  package  that  cannot  be  quantified  without 
a  tilt  sensor  (gyroscope  or  magnetometer).  These  uncorrected  motions  are  evident  in 
occasional  large  acceleration  spikes  that  cause  a  small  bias  (-'5-10  %)  in  wave  height 
estimates  (McIntyre,  2013). 

To  quantify  the  order  of  magnitude  of  the  errors  introduced  in  vertical 
acceleration  measurements  by  the  instrument  tilt  motions,  we  consider  the  idealized  case 
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of  an  accelerometer  package  rigidly  attaehed  to  a  surfaee-following  buoy.  As  the  buoy 
pitches  and  rolls  in  response  to  the  sea  surfaee  motion,  errors  are  introdueed  in  the 
vertical  acceleration  measurement  a^by  both  the  horizontal  aeeeleration  eomponents 

a^,  Qy  and  variations  in  the  statie  gravity  foree  eomponent  normal  to  the  sea  surfaee. 

These  latter  variations  have  the  same  time  scale  as  the  wave-indueed  surfaee  height 
fluctuations  and  thus  are  of  particular  concern  here.  For  small  values  of  the  surfaee  slopes 
and  the  error  (e)  in  the  measured  vertieal  aeeeleration  ean  be  expressed  as: 

£  =  a^ri^  +  -^(ril  +  ri^ ){a^  +  g)  (Al) 

In  the  linear  wave  theory  approximation  |<2j  =  g  and  at  the  surfaee 

’  «,v  =  -gVy  (A2a,  b) 

and  thus  Eq.  (Al)  reduees  to 

E  =  f  (');+');)  ■  (A3) 

This  analysis  shows  that  the  tilt  error  in  the  vertieal  aeeeleration  measurement  is 
eomparable  to  seeond  order  non-linear  effeets  and  thus  ean  be  formally  negleeted  in  the 
linear  approximation.  However,  the  quadratie  slope-squared  term  in  Eq.  (A3)  affeets  the 
measured  aeeelerations  not  at  the  primary  surfaee  wave  frequeneies,  but  at  the  sum  and 
differenee  frequencies  of  pairs  of  wave  eomponents.  Thus  in  mixed  swell-sea  wave 
fields,  that  are  eommonly  observed  in  the  oeean,  the  relatively  steep  wind  wave  slopes 
may  degrade  the  aeeeleration  measurements  at  lower  swell  frequeneies  where  is 

typically  much  smaller  than  g  .  This  sensitivity  is  eonfirmed  in  the  field  test,  whieh  show 
a  positive  bias  of  about  10%  in  swell  speetral  levels  from  the  aeeelerometer  data 
(Meintyre,  2013).  In  eontrast  to  the  expeeted  low  signal-to-noise  ratio  at  swell 
frequeneies,  the  measured  aeeelerations  of  steeper  wind  wave  eomponents  that  usually 
dominate  the  surfaee  slope  speetrum,  are  insensitive  to  this  souree  of  error.  Henee  the 
aeeelerometer  eomplements  the  GPS  sensor  that  is  better  suited  to  observing  the  longer 
period,  lower  steepness  swell  (Herbers  et  ah,  2012). 
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B,  VERTICAL  DISPLACEMENT  SPECTRA 

To  obtain  an  accurate  wave  spectrum  across  the  entire  swell  and  wind  wave 
bands,  we  form  a  eomposite  spectrum  using  the  GPS  estimate  below  0.15  Hz,  the 
aecelerometer  estimate  above  0.25  Hz,  and  a  weighted  average  in  between.  The  GPS 
(accelerometer)  weighting  was  varied  linearly  from  1(0)  at  0.15  Hz  to  0(1)  at  0.25  Hz. 
This  eomposite  wave  spectrum  produced  the  most  statistieally  reliable  estimates  of  H ^ 
when  eompared  to  the  Datawell  results.  The  resulting  WRD  speetra  are  in  excellent 
agreement  with  the  Datawell  speetra  (Figure  2).  Although  ground  truth  was  not  available 
for  direet  eomparison  with  the  speetral  tail,  the  speetrum  rolls  off  as  f~^ ,  whieh  is  in 
agreement  with  observations  of  equilibrium  wind-wave  spectra  (Toba  1973;  Phillips 
1985;  Battjes  et  al.  1987;  Alves  and  Banner  2003).  The  merged  speetral  estimates  of  the 
significant  wave  heights (with  E  the  variance  across  the  entire  0.05  -  1.0  Hz 
range)  were  found  to  be  less  than  3%  greater  than  those  of  the  Datawell. 

C.  SPATIAL  AVERAGING  ANALYSIS  TECHNIQUE 

The  deep-water  homogenous  dataset  was  analyzed  following  the  same  proeedures 
used  in  the  inhomogenous  tidal  inlet  region  in  order  to  demonstrate  the  aecuracy  of  the 
ensemble  averages  of  H ^  estimates  based  on  short  record  lengths  and  compared  with 
results  using  a  more  standard  20  minute  reeord  length.  As  expeeted,  the  varianee  in  H ^ 
values  estimated  using  the  same  short  segment  lengths  (410  s)  that  were  used  in  the  tidal 
inlet  analysis,  0.13  m  ,  is  larger  than  the  0.06  m  estimated  from  the  standard  20-minute 
segment  lengths.  However,  onee  the  spatial  and  ensemble  averaging  is  applied  to  the 
three  drifters,  the  variance  is  redueed  to  less  than  0.001  m  . 
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APPENDIX  B.  DIRECTIONAL  ANALYSIS  AND  BULK  STATISTICS 


A.  DIRECTIONAL  ANALYSIS 


The  direction  6  is  defined  in  a  coordinate  system  with  the  x  -axis  pointing  east 
and  the  y  -axis  pointing  north  ( 0  =0  for  waves  travelling  from  the  west,  90°  from  the 
south).  The  lowest  four  Fourier  moments  of  the  directional  distribution  of  wave  energy 
5(0)  at  each  frequency  can  be  expressed  in  terms  of  normalized  cross-spectra  (Long, 


1980) 
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where  C  and  Qare  the  co-  and  quadrature  part  of  the  cross-spectrum,  respectively. 


Equation  Bl  applies  to  x,y,z  displacements  but  can  be  adapted  to  velocity  and/or 


acceleration  measurements  by  using  transfer  functions  based  on  linear  wave  theory 
(Table  Bl). 


Estimates  of  0  and  dg  can  be  obtained  from  the  first-order  moments: 

tan(0)  =  —  (B2a) 

dg^^2(l-V?T^  (B2b) 

where  0  ,  defined  on  a  full  circle,  corresponds  to  the  direction  of  the  energy  flux  vector 
and  (Tg  can  be  interpreted  (in  the  limit  of  a  narrow  directional  spectrum)  as  the  standard 

deviation  of  5(0). 


Alternatively,  the  second-order  estimates  of  0  and  dg  are  given  by: 

tan(20)^^ 

«2 


(B3a) 
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(B3b) 

where  6  ,  defined  on  a  half  eirele,  eorresponds  to  the  major  axis  of  the  horizontal  wave 
orbital  displacements  and  dg  describes  the  associated  polarization  (e.g.  Herbers  et  ah, 
1999). 

Estimates  of  6  based  on  the  second-order  moments  (B3a)  are  generally  less 
useful  than  first-order  estimates  (B2a)  because  they  suffer  from  a  180°  ambiguity  in  the 
directional  estimates.  However,  the  second-order  6  and  dg  (B3a,  b)  depend  only  on  the 

horizontal  measurements  and  thus  provide  a  consistency  check  on  the  first-order 
estimates  that  use  the  vertical  acceleration  measurements  and  allow  for  comparing 
independent  estimates  based  on  GPS  position  and  velocity  data.  The  directional 
estimates  from  the  WRDs  are  in  fairly  good  agreement  with  the  Datawell  estimates 
(Figure  2).  Additionally,  the  first  order  moments  do  indeed  eliminate  the  180-degree 
ambiguity  that  is  evident  in  the  second-order  estimates  the  below  spectral  peak  and  near 
0.3  Hz. 

B,  BULK  WAVE  DIRECTIONAL  PARAMETERS 

In  addition  to  the  direction  spectral  estimates  6{f)  and  o[f),  estimates  of  bulk 

wave  parameters  were  obtained  from  short  wave  records  to  characterize  wave  conditions 
in  the  strongly  inhomogeneous  conditions  of  the  San  Francisco  Bight.  Energy-weighted 
bulk  estimates  of  6  and  o  were  obtained  by  integrating  the  auto-  and  cross-spectra  in 
Eq.  A1  over  the  swell-sea  frequency  range  (0.05-0.3  Hz).  The  bulk  directional  estimates 
were  found  to  be  within  5  degrees  of  the  Datawell  estimates. 
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TABLES 


SWAN 

WIND 

BOTTOM 

FRICTION 

CURRENTS 

Sim#l 

On 

On 

On 

Sim#2 

Off 

On 

On 

Sim#3 

On 

Off 

On 

Sim#4 

On 

On 

Off 

Table  1 .  Summary  of  SWAN  simulations  in  which  various  processes 

were  disabled. 


Displacement 

Spectra 

Eyy 

E. 

Q. 

Qy. 

WRD 

Euu 

Egg 

Egg 

C.g 

Egv 

(0^ 

0)^ 

(0^ 

-co^ 

-co^ 

-co^ 

Table  2.  Conversion  of  WRD  auto-  and  cross  spectra  of  vertical 
acceleration  [a)  and  horizontal  velocities  (m,v)  to  equivalent 

displacement  spectra  using  linear  theory  transfer  functions. 
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North  (UTM)  10 


FIGURES 


Figure  1.  Field  site  and  instrumentation.  The  field  site  is  in  the  San  Franeiseo 
Bight  where  30  drifters’  traeks  from  the  deployment  on  April  27,  2012, 
are  shown  in  blue.  A  rendering  of  the  WRD  and  a  photo  of  the  WRD 
floating  in  Monterey  Bay  are  shown. 
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Surface  Elevation  Spectrum  (m  /Hz) 


Figure  2.  Surface  height  spectral  estimates  from  measurements  made  in  the 
Monterey  Bay  by  a  Datawell  Buoy  (blue)  and  three  collocated  WRD 
buoys  (red).  Also  shown  is  the/“‘‘  spectral  roll  off  (black  dashed). 
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Figure  3.  Google  earth  image  of  the  SWAN  domain.  The  Point  Reyes  Buoy 
#46214  was  used  for  speetral  boundary  eonditions.  The  San  Franciseo 
Buoy  #46026  provided  wind  foreing  for  the  domain  outside  of  the  inlet. 
The  red  box  outlines  the  SWAN  domain.  The  small  white  box  indieates 
the  region  of  deereased  wind  foreing  (based  on  wind  and  water-level 

station  FTPCl). 
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Surface  Elevation  Spectrum  (m  /Hz) 


Figure  4.  Comparison  of  SWAN  spectrum  (red)  with  the  San  Francisco  Buoy 

spectrum  (blue)  at  1500  UTC. 
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Figure  5.  SWAN  significant  wave  height  and  mean  wave  direction  (a),  and  the 
DelftSD  depth  averaged  velocity  (b)  at  1515  UTC. 
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Figure  6. 
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Spatial  variability  of  significant  wave  height  from  the  April  27,  2012 
ebb  deployment.  Drifter  observations  (a)  are  compared  with  SWAN 
Sim#l  (b)  significant  wave  heights  that  are  evaluated  at  locations 
spatially  and  temporally  consistent  with  the  WRD.  The  WRD 
significant  wave  heights  were  estimated  using  the  merged  (GPS  and 
Accelerometer)  spectra.  The  estimates  were  obtained  from  6.83  min 
records  along  each  drift  track  and  averaged  spatially  over  a  500m  x  500 
m  grid.  The  regions  A-E  indicate  similar  wave  characteristics  used  for 
ensemble  averaging  regional  spectral  estimates  in  Figure  8a. 
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Relative  Frequency  (Hz) 


Figure  7.  April  27,  2012  Ensemble  averaged  surface  height  spectra  in  regions 
A-E  (A=offshore,  B=over  the  SE  Bar,  C=east  of  SF  Bar,  D=approaches 
to  SF  Inlet,  E^inlet,  See  Figure  6).  Observed  spectra  (top)  are  compared 
with  SWAN  Sim  #1  spectra  calculated  at  the  mean  time  of  the  drifter 
ensemble  (bottom).  Also  indicated  (top)  are  the  f~'^  spectral  roll  off 
(black  dashed)  and  the  blocking  frequency  in  regions  B-E. 
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Figure  8.  WRD  significant  wave  height  estimates  (color)  with  ray- trajectories 
of  1 1  second  swell  (black  lines)  on  bathymetry  (grey  curves).  Rays 
were  calculated  with  wave  current  interaction  using  depth-averaged 
currents  from  DelftSD  (a),  and  with  no  currents  (b). 


54 


4.19 


H  (m) 


4.17^ - ^ 

5.31  5.38  5.45 


(£> 

o 


o 

z 


4.19 


4.18 


4.17 


I 


5.31  5.38  5.45 

East  (UTM)  10® 


3 

2.5 
2 

1.5 
1 


Figure  9.  Significant  wave  heights  predicted  in  SWAN  Sim  #2  in  which  wind 
is  disabled  (a),  SWAN  Sim#3  with  bottom  friction  disabled  (b),  and 
SWAN  Sim  #4  with  currents  disabled  (c). 
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Figure  10. 
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Ray  trajectories  of  1 1  second  swell  (accounting  for  depth  and  current 
refraction)  superposed  on  WRD  estimates  of  bulk  directional  spread 
(second-order  estimates,  Appendix  B). 
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Figure  1 1 .  WRD  significant  wave  height  estimates  (color)  with  ray- trajectories 
of  5  second  wind  wave  (black  lines)  on  bathymetry  (grey  curves).  Rays 
were  calculated  with  wave  current  interaction  using  depth-averaged 
currents  from  DelftSD  (top),  and  with  no  currents  (bottom). 
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Figure  12.  Directional  spectral  estimates  from  the  Datawell  buoy  (blue)  and  the 
WRD  (red)  deployed  alongside  the  Datawell  buoy  during  the  Monterey 
Bay  experiment.  From  top  to  bottom;  mean  direction  and  directional 
spread  (as  functions  of  frequency),  respectively.  The  left  panels  show 
“traditional”  estimates  based  on  measurements  of  both  vertical  and 
horizontal  surface  motions.  The  right  panels  show  alternative  estimates 
based  solely  on  horizontal  surface  motions. 
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Figure  13.  Google  Earth  image  of  the  Mouth  of  the  Columbia  River  and 

offshore  approaches.  The  yellow  box  indicates  the  SWAN  domain  and 
the  red  box  outlines  the  jetties  at  the  entrance  to  the  Columbia  River. 
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Figure  14.  Polar  spectral  density  plot  from  CDIP  station  179  in  Astoria  Canyon. 

This  location  is  representative  of  the  wave  conditions  offshore  of  the 
Mouth  of  the  Columbia  River,  http://cdip.ucsd.edu 
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Figure  15.  Three  different  types  of  Wave  Resolving  Drifters  (WRD).  The 

WRD-A  deployed  from  the  R/V  QUESTUARY  (top  left).  The  mix  of 
WRD  A-B-J  on  the  deck  of  the  R/V  POINT  SUR  (top  center).  The 
WRD  boat  hook  recovery  (top  right).  The  more  advanced  WRD-J 

(bottom). 
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Figure  16.  Drift  velocity  measurements  along  drift  tracks  averaged  temporally 
over  204.8  seconds  and  spatially  onto  a  100m  x  100m  grid.  The 
measurements  are  plotted  on  the  centroid  of  the  data  within  each  grid 
cell.  The  color  bar  on  the  right  indicates  the  current  scale.  The 
rectangular  boxes  indicate  regions  where  wave  spectra  were  evaluated. 
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Figure  17.  Spatial  variability  of  significant  wave  height  from  the  June  8,  2013 
ebb  deployment.  The  WRD  observed  significant  wave  heights  were 
estimated  using  the  combined  GPS-Accelerometer  spectra.  The 
estimates  were  obtained  from  3.4  min  records  along  each  drift  track  and 
averaged  spatially  over  a  100m  x  50  m  grid.  The  colored  zones  coincide 
with  the  regional  spectral  estimates  in  Figure  19. 
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Figure  18.  June  08,  2013  Ensemble  averaged  surface  height  spectra  in  Zones 
03-13  (03=offshore,  02=over  the  mound,  01=east  of  mound,  B=over 
the  Bar,  ll=inside  the  Bar,  12=in  the  inlet,  13=in  the  river.)  Surface 
elevation  spectrum  (a),  mean  direction  (b),  and  directional  spread  (c). 
The  directional  estimates  are  based  on  2"‘*  order  moments  (Appendix 

B). 
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Figure  19.  Regional  histograms  of  the  individual  WRD  significant  wave  height 
estimates  within  each  zone.  The  zones  are  color  coded  to  match  the 
spatial  regions  and  spectra  in  Figure  18. 
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Figure  20.  Zones  of  interest  (eolor)  with  ray- trajectories  of  14-second  swell 
waves  from  an  offshore  angle  of  280  degrees  (black  lines)  on 
bathymetry  (grey  curves).  Rays  were  calculated  with  wave  current 
interaction  using  depth-averaged  currents  from  SELFE. 
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Figure  21 .  Zones  of  interest  (eolor)  with  ray- trajectories  of  14-second  swell 
waves  from  an  offshore  angle  of  270  degrees  (black  lines)  on 
bathymetry  (grey  curves).  Rays  were  calculated  with  wave  current 
interaction  using  depth-averaged  currents  from  SELFE. 


67 


North  (UTM)  10 


5.124 


5.12 


5.116 


4.04 


4.08 


4.16 


East  (UTM)  10' 


Figure  22.  Zones  of  interest  (color)  with  ray-trajectories  of  5 -second  wind 
waves  from  an  offshore  angle  of  285  degrees  (black  lines)  on 
bathymetry  (grey  curves).  Rays  were  calculated  with  wave  current 
interaction  using  depth-averaged  currents  from  SELFE. 
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Figure  23.  Ray-trajectories  of  5 -second  wind  waves  from  an  offshore  angle  of 

285  degrees  (black  lines)  on  bathymetry  (grey  curves).  Rays  were 
calculated  with  wave  current  interaction  disabled. 
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Figure  24.  SWAN  Sim  #1  significant  wave  height  predictions.  The  Zones  of 
interest  coincide  with  those  from  the  WRD  discussion  and  the  SWAN 

spectra  in  Figure  25. 
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Figure  25.  SWAN  spectra  from  Sim  #1  are  calculated  at  grid  points  near  the 

center  of  the  Zones. 
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Figure  26.  Comparison  of  significant  wave  height  estimates  from  the  WRD 

estimates  and  SWAN  Sim  #1  (a)  and  Sim  #4  (b)  collocated  in  space  and 
time.  The  colors  are  consistent  with  the  Zones  in  Figures  24  and  25. 
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Figure  27.  SWAN  significant  wave  height  with  currents  disabled. 
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Figure  28.  Significant  wave  height  estimates  from  the  WRD  with  rays  from  an 
offshore  wave  direction  of  275  and  a  14  second  period. 
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Figure  29.  Ensemble  averaged  surface  height  spectra  in  Zone  B  for  enhanced 
(red)  and  un-enhanced  (blue)  local  wave  conditions.  Surface  elevation 
spectrum  (a),  mean  direction  (b),  and  directional  spread  (c). 
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Figure  30.  Significant  wave  height  plotted  against  directional  spread  in  Zone  B. 
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tgure  31.  A  time  series  of  orbital  veloeities  from  one  of  the  drifters  pa: 
over  the  Bar.  The  time  series  is  shown  in  reverse  to  match  the  s| 

variation  in  Figure  29  from  4. 12  to  4.14  10^  (UTM). 
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